A REVIEW  OF  RESISTANC  rlU  IW*-?  .pOo  JNL 


<2^ 


DAVID  W.  TAYLDR  NAVAL  SHIP 
RESEARCH  AND  DEVELOPMENT  CENTER 


1 

3 


B«thMda.Md.  20064 


A REVIEW  OF  RESISTANCE  PREDICTION 
AND  DESIGN  METHODS  FOR  STRUT-POD  INLETS 


Pierre  Lafrance 


^ APPROVED  FOR  PUBLIC  RELEASE:  DISTRIBUTION  UNLIMITED 


SHIP  PERFORMANCE  DEPARTMENT 
DEPARTMENTAL  REPORT 


March  1978 


SPD-  789-01 


! 


MAJOR  DTNSRDC  ORGANIZATIONAL  COMPONENTS 


■■VI 


;gCuniTV  CLASSIFICATION  OF  THIS  FAOC  (Whan  Otm  BnlwMQ 


REPORT  DOCUMENTATION  PAGE 


|2.  OOVT  ACCESSION  NO. 


789-61 


^ JEVIEW  OF  BESISTANCE^REDICTION  AND  DESIGN 
^ETHODS  FOR  JTRUT-^OD  JnLETS^  ^ ^ 


IT.  AUTHONfaJ 


Pi erre/Laf ranee 


tr  FCNFONMINO  ONOANIZATION  NAME  ANO  AODNESS 

David  W.  Taylor  Naval  Ship  R&O  Center 
Bethesda , MO  20084 

II.  CONTNOLLINO  OFFICE  NAME  ANO  AOONESS 


READ  nrafftOCTKNIS 
BEFORE  OOllPLBTIWO  rOKM 

S.  NECIFIENT'S  CATALOG  NUMEEN 


S.  TVFE  OF  REFONT  * PERIOD  COVERED 


t.  RERFDRMINO  ORG.  REPORT  NUMBER 


CONTRACT  OR  grant  NUMBBRfaJ 


Naval  Sea  Systems  Command 
Washington,  D.C.  20350 


10.  PROGRAM  ELEMENT.  PPOJECT,  TASK 
AREA  A WORK  UNIT  NUMBERS 

Element  62501.^ 

Task  12501  1 

Task  Area/S^34323il  1 

iz.  report’otWI^ * 

/ Mar  A W78  ^ 

IS.  NUMBER  Of  PAGES 


monitoring  agency  name  * AOORESSrir  dUi0rm»t  trom  Contt^lUng  Otftc*)  IS.  SECURITY  CLASS,  (oi  thl9  nporl) 


UNCUSSIFIED 


11a.  oeclassification/downgrading 
SCHEDULE 


IS.  distribution  statement  (ot  thlt  RapofU 


APPROVED  FOR  PUBLIC  RELEASE:  DISTRIBUTION  UNLIMITED 


I IT.  DISTRIBUTION  STATEMENT  (ot  Ihm  obotemet  ontotod  In  Block  20,  It  ditloront  trom  Rai 


IS.  SURRLEMEnTARY  NOTES 


I If.  KEY  WORDS  (Continuo  on  roootoo  oldo  It  noeooomrf  m4  Idontlty  by  block  ntmbor) 


Waterjet,  Strut-Pod,  Inlet,  Drag 


20.  abstract  (ConihHto  on  roroeoo  oldo  1/  nocoooory  md  Idonitfy  by  block  mmbot) 

'^Two  aspects  of  strut-pod  waterjet  inlet  technology  are  reviewed  in 
this  report.  First,  hydrodynamic  aspects  of  inlet  design  are  discussed 
and  various  design  methods  are  presented.  Second,  the  subject  of  drag 
prediction  is  briefly  reviewed,  state-of-the-art  drag  prediction  methods 
are  presented  and  their  predictive  accuracy  is  compared. 

For  superhigh  speed  applications  (0.^100  kts),  a variable  area  concept 
with  auxiliary  strut  vanes  is  recommended.  For  speeds  below  50  kts,  fixed 


DO  1473  EDITION  OF  I NOV  EE  IE  OBEOLETk  ^ 

$/N  0103.LF414.4E0I 


’ ^wmssiFiEP  tn , 

SECURITY  CLAEEIFICATION  OF  THIS  PAGE  rWAan  Oatn 


UNCLASSIFIED 


»tCUWTV  CL^MIFICATIOM  or  this  HAOI  rmtm  0«a  BntwMI) 


20.  Continued 


area  Inlets  should  suffice.  Intermediate  speed  applications  should  be 
Investigated  on  a case-by>case  basis  to  determine  If  auxiliary  vanes 
are  desirable.  . 


UNCLASSIFIED 


TABLE  OF  CONTENTS 


ABSTRACT 


ADMINISTRATIVE  INFORMATION 


INTRODUCTION 


DESIGN  PHILOSOPHY 


PRACTICAL  DESIGN  CONSIDERATIONS 


DESCRIPTION  OF  EXISTING  DESIGN  METHODS 


RAM  INLET  PERFORMANCE 


DISCUSSION  AND  RECOMMENDATION  FOR  INLET  DESIGN 


ANALYTICAL  DRAG  PREDICTIONS 


SKIN  FRICTION  DRAG 


CAVITY  DRAG 


INTERFERENCE  DRAG 


PREDICTIVE  STATE-OF-THE-ART 

DRAG  PREDICTION  METHODS  ACTUALLY  USED  IN  DESIGN  PROGRAMS 


REVIEW  OF  THE  EXPERIMENTAL  DATA 


CONCLUSION  AND  RECOMMENDATION  CONCERNING  INLET  DRAG 


RECOMMENDATIONS  FOR  FURTHER  WORK 


REFERENCES 


Section  B' 
g.;V  Section  • D 


ACKNOWLEDGEMENT 


LIST  OF  FIGURES 


Page 


Figure  1 - Waterjet  Propulsion  Unit 4 

Figure  2 - Pod  with  Slat  and  Auxiliary  Doors 5 

Figure  3 - Flow  In  the  Vicinity  of  the  Inlet 8 

Figure  4 - Variation  of  Inflow  Angle  with  Inlet  Velocity  Ratio 

(IVR)  for  a Typical  Inlet 9 

Figure  5A  - Effect  of  Inlet  Loss  Coefficient  on  Waterjet 

Optimization 15 

Figure  5B  - Effect  of  Inlet  Drag  Coefficient  on  Waterjet 

Optimization.. 16 

Figure  6A  - Variation  in  Operational  Inlet  Velocity  Ratio  with 
Vehicle  Speed  at  a Constant  Flow  Rate,  Design 
Speed  = 5D  kt 19 

Figure  6B  - Variation  In  Operational  Inlet  Velocity  Ratio  with 
Vehicle  Speed  at  a Constant  Flow  Rate,  Design 
Speed  = 100  kt 20 

Figure  7 - Variation  in  Physical  Size  and  Contour  of  Subcavitating 

Inlet  Lips 21 

Figure  8 - Variation  in  Physical  Size  of  Subcavitating  Inlet 

Nacelles  Designed  for  a Speed  of  100  Knots  and  a 
Constant  Flow  Rate  of  23,500  Ib/sec  (10660  kg/sec)  with 
Design  Inlet  Velocity  Ratio 22 

Figure  9 - Cavitation  Characteristics  of  the  Six  Inlet  Designs 24 

Figure  lOA-  Variation  with  Design  Inlet  Velocity  Ratio  of  the 
Components  of  the  Total  External  Drag  for  the 
Subcavitating  Inlet  System 25 

Figure  lOB-  Variation  with  Design  Inlet  Velocity  Ratio  of  the 

Components  of  the  Total  External  Drag  of  the  Base  Vented 
(at  Maximum  Thickness)  Inlet  System 26 

Figure  IOC-  Variation  with  Design  Inlet  Velocity  Ratio  of  the 

Components  of  the  Total  External  Drag  of  the  Base  Vented 
(Forward  of  Maximum  Thickness)  Inlet  System 27 

Figure  lOD-  Variation  with  Design  Inlet  Velocity  Ratio  of  the 
Components  of  the  Total  External  Drag  of  the 
Supercavltating  System 28 

Figure  11  - Comparison  of  the  Total  External  Drag  Coefficients  ot 
the  Inlet  Systems  Investigated  vs.  Design  Inlet 
Velocity  Ratio 29 


■ I 

r 


1v 


Page 


[ 


> 

[ 


r 


j 

* 

E 

\ 


I 

1 

I, 

I 


f 


Figure  12  - Total  Inlet  Loss  Coefficient  vs.  Design  Inlet 

Velocity  Ratio 30 

Figure  13  - Component  Head  Loss  vs.  Flowrate  for  a Design  Inlet 

Velocity  Ratio  = 0.75 31 

Figure  14A-  Optimum  Design  Inlet  Velocity  Ratio  for  Subcavitating 

Inlet 33 

Figure  14B-  Optimum  Design  Inlet  Velocity  Ratio  for  Inlet  Base 

Vented  at  the  Point  of  Maximum  Thickness 34 

Figure  14C-  Optimum  Design  Inlet  Velocity  Ratio  for  the  Inlet  Base 

Vented  Forward  of  the  Point  of  Maximum  Thickness 35 

Figure  14D-  Optimum  Design  Inlet  Velocity  Ratio  for  Supercavitating 

Inlet 36 


Figure  15  - Configuration  in  the  Retracted  Position  of  Center  Body..  40 


Figure  16  - 
Figure  17  - 
Figure  18  - 

Figure  19  - 
Figure  20  - 


Internal  Pressure  Loss  Performance  - High  Speed 
Conditions 41 

Internal  Pressure  Loss  Performance  - Low  Speed 
Conditions 42 

Inlet  Cavitation  Boundaries  Measured  and  Predicted  for 
the  34.73  Inch  (0.882  m)  Centerbody  Extension  (Fully 
Extended) 44 


Inlet  Cavitation  Boundaries,  Measured  and  Predicted, 

for  the  0.0  Inch  Centerbody  Extension  (Fully  Retracted).  45 

Typical  Inlet  Drag  Performance  for  a Range  of  IVR, 

Centerbody  Retracted,  Simulated  20  Knots 47 


Figure  21  - 

Figure  22  - 
Figure  23  - 
Figure  24  - 
Figure  25  - 
Figure  26  - 


Typical  Inlet  Drag  Performance  for  a Range  of  IVR, 

Centerbody  Extended,  Simulated  100  Knots 48 

Subcavitating  Model  - Overall  Pressure  Loss 54 

Drag  Coefficient  for  Subcavitating  Lockheed  Model 55 

Supercavitating  Model  - Overall  Pressure  Loss 57 

Drag  Coefficient  for  the  Supercavitating  Lockheed  Model.  58 

Comparison  of  Cavitation  Characteristics  for  Primary 
Inlet  to  Pod  Diameter  Ratios  of  40%,  45%,  and  50% 60 


L 


V 


Page 


I 


Figure  27  - Internal  Cavitation  Inception  Speed  vs.  Effective 
Inlet  Velocity  Ratio  for  45%  Inlet  Covrl;  Constant 
Auxiliary  Inlet  Opening  of  8%  Inches  (0.216  tn) 62 

Figure  28  - Sunmary  Comparison  for  Computed  and  Experimental 
Residuary  Drag  Coefficients  (C-  = Cr.  - C,)  for  all 
Cases 80 

Figure  29  - Nondimenslonal  Plot  (x/L,  y/L)  of  Series  4620  - 2, 

3,  4 Models 82 


LIST  OF  TABLES 


TABLE  1 - RESULTS  OF  DESIGN  INLET  VELOCITY  OPTIMIZATION 37  , 

TABLE  2 - SU^f1ARY  OF  INLET  PERFORMANCE 64 

TABLE  3 - SKIN  FRICTION  COEFFICIENTS 77 

TABLE  4 - RMS  % ERROR  IN  HESS'  DRAG  COMPUTATIONS 78 

TABLE  5 - APPROXIMATE  PREDICTIVE  ACCURACY  OF  DRAG  PREDICTION 

METHODS 100 

I 


I 

I 


Symbol 


Meaninc 


Area,  Inlet  area 


NOMENCLATURE 

Dimension 

.2 


Chord 

Drag  coefficient 
Skin  friction  coefficient 
Pressure  coefficient 
Residuary  drag  coefficient 
Inlet  diameter 


Straight  pipe  friction  factor  - 
Fanning  friction  factor 
Gravitational  acceleration  LT' 


Boundary  layer  shape  factor 
Inlet  head  losses 
Submergence  depth 
Inlet  velocity  ratio  = V^/V^ 
Inlet  loss  coefficient 
Total  head  loss  coefficient 


l>0.  Pi 


Body  length  L 

Reynolds  number 

-1  -2 

Pressure,  inlet  pressure  ML  T 

-1  -2 

Dynamic  pressure  ML  T 

-1  -2 

Free  stream  stagnation  pressure  ML  T 


Volumetric  flowrate 


l3t-1 


Radial  coordinate 


vii 


Meaning 

Dimension 

Radial  coordinate  of  a point  on  body 

L 

Strut  thickness 

L 

Free  stream  velocity 

LT' 

Streamwise  velocity 

LT‘ 

Surface  velocity 

LT' 

Fluid  speed 

LT' 

Inlet  flow  speed 

LT' 

Jet  flow  speed 

LT' 

Craft  speed 

LT' 

Coordinate  measured  along  body 

L 

Coordinate  measured  normal  to  body 

L 

Diffuser  expansion  half-angle 

- 

Boundary  layer  thickness,  where  u = 0.99  u^ 

00 

L 

Boundary  layer  displacement  thickness  = (1  - 

.00 

iL.)ILdy 

00  0 

L 

Boundary  layer  momentum  thickness  = (1  - 

*'0 

y_ILdy 

00  0 

L 

Jet  propulsive  efficiency 

- 

Fluid  kinematic  viscosity 

L^T' 

Fluid  density 

ML'' 

Cavitation  number 

- 

Stagnation  streamline 

- 

I [ Two  aspects  of  strut-pod  waterjet  inlet  technology  are  reviewed  in 

! I this  report.  First,  hydrodynamic  aspects  of  inlet  design  are  discussed 

I I and  various  design  methods  are  presented.  Second,  the  subject  of  drag 

. I prediction  is  briefly  reviewed,  state-of-the-art  drag  prediction  methods 

: [ are  presented  and  their  predictive  accuracy  is  compared. 

; For  superhigh  speed  applications  (=  100  kts),  a variable  area 

I concept  with  auxiliary  strut  vanes  is  recommended.  For  speeds  below  50 

I 1 kts,  fixed  area  inlets  should  suffice.  Intermediate  speed  applications 

i f should  be  investigated  on  a case-by-case  basis  to  determine  if  auxiliary 

i I vanes  are  desirable. 

i I 

f I 

1 ADMINISTRATIVE  INFORMATION 


This  project  was  sponsored  by  the  Naval  Sea  Systems  Command,  Code 
0331G,  Task  Area  SF  43432301,  Task  12501,  Element  62543N. 


INTRODUCTION 


The  design  of  inlets  for  waterjet-propelled  craft  shares  a common 
feature  with  many  other  technology-intensive  fields.  It  has  reached  the 
stage  where  most  of  the  knowledge  is  scattered  throughout  a multi tute  of 
reports  and  technical  papers  from  which  a comprehensive  view  is  difficult 
to  obtain.  It  would  appear  that  a compendium  of  inlet  design  methodology 
would  be  a timely  contribution.  However,  the  very  speed  with  which 
developments  are  being  made  in  engineering  and  scientific  fields  touching 
on  that  area  precludes  the  present  realization  of  a work  of  enduring 
value.  However,  it  is  safe  to  state  that  the  major  problems  in  inlet 
design,  at  least  in  the  hydrodynamic  sense,  have  been  Identified,  and 


that  a document  describing  these  problems  and  the  various  existing 
solutions  is  not  out  of  place. 


This  report  represents  a step  toward  the  establishment  of  an  inlet 
design  philosophy  at  DTNSRDC.  Although  a detailed  mathematical 
description  of  the  hydrodynamics  pertaining  to  inlet  design  falls  outside 
the  scope  of  this  report,  the  underlying  physical  principles  are 
presented  to  give  the  reader  an  appreciation  for  the  design  methods  which 
have  evolved  over  the  years. 

The  actual  design  of  ram  inlets  is  reviewed  through  the 
descriptions  of  several  existing  design  methods.  These  descriptions  are 
necessarily  brief  and  are  meant  to  be  a short  review  of  the 
state-of-the-art.  The  reader  is  led  to  the  references  for  additional 
details. 

Since  the  purpose  of  this  study  was  to  gather  from  the  literature 
examples  representative  of  the  current  state-of-the-art,  some  passages  in 
this  report  have  been  transcribed  almost  verbatim  from  existing  reports 
and  papers.  These  are  not  identified  as  direct  quotations,  but  the 
sources  are  all  listed  in  the  references. 

DESIGN  PHILOSOPHY 

In  principle,  there  is  little  difference  between  conventional  propellers 
and  waterjet  propulsion  systems.  With  the  former,  the  actuator  is 
situated  in  the  external  stream  while  with  the  latter,  the  stream  is 
directed  within  the  hull,  where  it  is  allowed  to  interact  with  the 
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actuator  and  is  subsequently  discharged  astern  in  the  form  of  a high 
speed  jet.  There,  the  shape  of  the  ducting  dictates  the  flow  speed  and 
hence  can  be  advantageously  used  to  increase  static  pressure  at  the  pump, 
thereby  retarding  cavitation.  Waterjet  inlet  design  centers  around  the 
determination  of  ducting  configurations  which  will  adequately  supply  the 
pump  while  maintaining  hydrodynamic  losses  at  a minimum.  The  principal 
difference,  from  a designer's  point  of  view,  between  propellers  and 
waterjet  propulsion  is  that  with  waterjets, one  has  more  control  over  the 
incoming  flow. 

A waterjet  propulsion  system  is  built  around  a central  element  - 
the  pump  - whose  purpose  is  to  increase  the  momentum  of  the  water  stream 
flowing  through  it,  resulting  in  a net  thrust  on  the  craft.  A concise 
and  lucid  description  of  waterjet  pumps  was  given  by  Wislicenus^T  Figure 
1 schematically  shows  a typical  waterjet  propulsion  unit  utilizing  a pod 
inlet.  Pump  performance  is  limited  by  cavitation.  To  operate  properly, 
the  pump  must  be  supplied  with  the  required  flowrate  at  pressures  which 
will  not  result  in  cavitation.  In  addition,  the  flow  at  the  pump  face 
should  be  uniform.  Providing  the  pump  with  these  flow  conditions  is 
accomplished  by  the  conduit  upstream  of  the  pump  face.  This  structure 
typically  comprises  the  inlet,  diffuser,  turning  vanes,  ducting,  and  may 
include  devices  for  varying  the  effective  inlet  area,  such  as  auxiliary 
inlets  or  a retractable  centerbody.  These  elements  are  housed  within  the 
nacelle  and  strut,  as  shown  in  Figure  2,  and  shall  be  referred  to  by  the 

generic  names  of  "waterjet  inlet,"  or  simply  "inlet." 

s 

1 

♦List  of  references  given  on  page  104  j 
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The  purpose  of  the  Inlet  1s  to  match  the  external  stream  to  the 


pump,  maintaining  cavitation  - and  separation  - free  operation  at 

specified  operating  conditions  (e.g.,  hump  and  cruise  speeds)  while  • 

keeping  losses  and  system  weight  at  a minimum.  Since  large  losses  are 

likely  to  occur  In  the  Inlet,  It  Is  Important  that  this  area  be  designed 

with  great  care. 

In  general,  for  a hydrofoil  vehicle  a circular  ram  Inlet  located  at 

2 

the  strut-foil  Intersection  Is  believed  to  offer  the  greatest  advantage  . 

The  ram  Inlet,  considered  a very  efficient  Inlet  configuration,  achieves 
a high  pressure  recovery  and  low  drag.  The  Inlet  total  pressure  may  be 
In  the  order  of  98%  of  free  stream  total  pressure.  In  addition,  the 
circular  Inlet  and  nacelle  offer  a minimum  circumference  and  therefore 
the  least  surface  drag  for  a given  Inlet  area.  * 

Waterjet  systems,  to  be  fair,  are  not  free  of  disadvantages.  A 
substantial  price  Is  paid  In  inlet  losses  and  added  drag.  In  addition, 
total  Installed  weights  of  waterjet  systems  tend  to  be  greater  than 
similar  propeller  Installations.  The  options  and  tradeoffs  to  be 

considered  In  choosing  between  these  two  propulsion  systems  for  a < 

particular  application  constitute  an  Interesting  subject  which  lies 
outside  the  scope  of  this  study. 

* It  should  be  noted  that  rectangular  inlets  can  be  used  to  simplify 
the  transition  to  the  strut  duct  and  thereby  minimize  internal  losses 


PRACTICAL  DESIGN  CONSIDERATIONS 


« The  design  of  a waterjet  Inlet  presents  a formidable  hydrodynamics 

problem  touching  upon  several  areas  where  fundamental  understanding  Is 
still  lacking.  In  response  to  the  Immediate  need  for  practical  designs, 
several  approximations  have  been  developed  to  simplify  otherwise 
cumbersome  procedures  as  well  as  to  empirically  compensate  for  a lack  of 
understanding  in  certain  areas.  The  art  of  inlet  design  has  evolved 
around  these  approximations  and  the  ways  in  which  they  are  incorporated 
Into  various  design  programs. 

In  this  section,  we  shall  examine  the  various  components  of 
waterjet  pod  inlets  and  discuss  the  factors  which  must  be  considered  in 
their  design.  This  discussion  will  be  maintained  on  a very  general 
level,  the  details  to  be  saved  for  the  following  section,  where  we 

i 

I 

j • discuss  particular  design  procedures. 

i For  waterjet  inlets,  and  especially  for  inlets  housed  in  a pod,  we 

I ' 

must  distinguish  between  two  separate  flow  regions.  By  internal  flow,  we 
shall  denote  the  flow  bounded  by  the  centerline  and  the  stagnation 
streamline  as  shown  in  Figure  3.  The  internal  flow  passes  through 
the  pump  and  is  ejected  as  a high  speed  jet.  The  external  flow  is  also 
bounded  by  and  passes  over  the  external  surface  of  the  pod. 

The  basic  criteria  governing  the  design  of  external  nacelle  shape 
and  strut  are  low  drag  and  cavitation- free  performance.  This  places 
limitations  on  total  wetted  surface  area  and  on  the  maximum  surface  ) 

velocity.  At  very  high  craft  speed  (say  100  knots),  it  may  become 
advantageous  to  ventilate  parts  of  the  strut-pod  combination  to  reduce 
t ' the  wetted  area. 
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Figure  3 - Flow  in  the  Vicinity  of  the  Inlet 

The  inlet  nacelle  is  assumed  to  be  sufficiently  removed  from  the 
hull  that  the  latter's  boundary  layer  has  negligible  effect  on  the 
former.  The  inlet's  effect  on  hull  pressure  distribution  can  affect 
overall  craft  performance.  Conversely,  the  hull,  to  a degree  dependent 
upon  strut  length,  influences  the  flow  in  the  vicinity  of  the  inlet. 

These  effects  are  usually  not  considered  in  the  design  of  strut-pod 
Inlets. 

The  inlet  lip  is  a very  critical  area  of  the  inlet  as  it  is 
extremely  susceptible  to  cavitation  and  contributes  significantly  to 
drag.  Figure  4 shows  the  flow  in  the  vicinity  of  the  lip  for  three 
values  of  the  inlet  velocity  ratio.  At  low  IVR's,  the  stagnation  point 
is  located  on  the  interior  side  of  the  lip,  forcing  the  external  flow  to 
turn  around  the  lip,  where  it  may  experience  large  accelerations 
depending  on  the  severity  of  the  leading  edge  curvature.  If  sufficiently 
high  flow  speeds  are  reached,  the  pressure  will  drop  to  very  low  values 
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Figure  4 - Variation  of  Inflow  Angle  With  Inlet  Velocity  Ratio 
(IVR)  for  a Typical  Inlet 
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and  cavitation  will  appear  on  the  external  part  of  the  lip.  At  high 
IVR's  the  stagnation  point  occurs  on  the  external  portion  of  the  lip  and 
consequently,  cavitation  Is  likely  to  manifest  Itself  on  the  Inside  lip 
surface.  This  can  have  a particularly  deleterious  effect  on  performance 
since  with  Increasing  severity.  Internal  cavitation  can  choke  the  Inlet 
and  starve  the  pump.  The  obvious  answersto  lip  cavitation  are,  (1)  to 
operate  near  IVR's  as  near  unity  as  possible,  (2)  to  minimize  the  flow 
accelerations  by  Increaslnq  the  lip  leading  edge  radius,  and  (3)  variable 
geometry.  These  remedies  are  not  without  difficulties,  however.  Keeping 
the  IVR  near  unity  over  a broad  range  of  operating  speeds  may  necessitate 
the  use  of  variable  geometry  while  decreasing  leading  edge  curvature  is 
accompanied  by  an  Increase  In  drag. 

While  leading  edge  radius  Is  a fixed  quantity  In  a given  design, 
several  types  of  variable  Inlet  geometry  have  been  considered  to 
accomnodate  greater  flowrates  at  hump  speed  without  risking  Inlet 
cavitation,  two  of  which  are  shown  In  Figure  2. 

The  IVR  Is  an  Important  parameter  1n  Inlet  design.  An  inlet 
velocity  ratio  of  unity  should  be  selected  for  the  design  cruise 
condition.  From  practical  considerations  however,  (for  a fixed  area 
inlet)  this  value  should  usually  be  somewhat  less  than  1.0.  For  example, 
the  speed  of  a hydrofoil  craft  at  hump  is  approximately  one  half  the 
cruise  speed.  For  the  same  flowrate,  unit  IVR  at  cruise  implies  an  IVR 
of  2 at  hump,  making  cavitation  on  the  Internal  lip  very  likely. 
Accordingly,  the  cruise  IVR  must  be  lowered  as  much  as  possible  without 
incurring  cavitation  on  the  external  lip.  A viable  solution  to  this 
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design  problem  Is  the  use  of  auxiliary  Inlets  or  variable  geometry  to 
provide  additional  Inlet  area  at  hump  with  a reduction  In  Inlet  velocity. 

■ V Merging  with  the  Inlet  lip  and  extending  to  the  pump  face  Is  a duct 

section  called  the  diffuser.  For  the  purpose  of  analysis,  the  diffuser 
I may  be  viewed  as  comprising  three  separate  components,  the  Inlet 

r diffuser,  the  elbow,  and  the  strut  diffuser.  The  primary  function  of  the 

I diffuser  Is  to  supply  the  pump  with  the  required  flowrate  while 

^ maintaining  a minimum  total  pressure  loss  and  flow  distortion. 

I Calculation  of  energy  losses  in  such  a complex  system  is  well  beyond 

present  analytical  capabilities  and  consequently,  heavy  use  is  made  of 
empirical  results.  Diffuser  performance  is  characteristically  expressed 
in  a "diffuser  recovery  factor"  or  pressure  recovery,  which  is  a measure 
of  the  recoverable  ingested  kinetic  energy.  The  remaining  energy 
fraction  is  irreversibly  dissipated  by  such  processes  as  skin  friction  or 
• the  generation  of  turbulent  eddies  following  boundary  layer  separation, 

and  contributes  to  inlet  losses. 

Diffuser  contour  shapes  greatly  affect  performance  and  the  bell 
shape  contour  has  been  shown  to  produce  higher  pressure  recovery  than 

I 3 

equivalent  straight  wall  or  trumpet  shaped  diffusers.  The  reason  given 

I 

' for  this  is  that  the  largest  positive  pressure  gradient  occurs  where  the 

i 

boundary  layer  Is  still  thin  and  able  to  overcome  the  retarding 
influence. 


A study  of  prerotation  of  the  flow  just  ahead  of  the  pump  to 
Increase  suction  performance  and  ensure  maximum  takeoff  thrust  at  low 
craft  speeds  has  recently  been  experimentally  Investigated  by  Rockwell 


4 

International.  The  purpose  of  this  study  was  to  Investigate  waterjet 
suction  performance  gains  obtainable  with  a system  of  prerotator  vanes. 

A set  of  vanes,  providing  the  requisite  solid  body  prerotation,  was 
designed,  fabricated  and  tested  with  an  existing  scale-model  Powerjet  20 
Inducer  In  a water  tunnel.  Analysis  of  test  results  Indicate  that 
Introducing  prerotation  by  means  of  an  upstream  vane  system  may  likely 
create  losses  which  will  detract  from  gain  1n  inducer  suction  performance 
obtained  purely  from  flow  prerotation.  From  pressure  measurements, 
upstream  pressure  losses  were  found  to  be  significant  and  possibly  caused 
by  vane  blockage  and  flow  separation.  Further  tests  were  recommended 
with  a new  vane  geometry. 

The  problem  of  avoiding  cavitation  in  the  duct  system  is  most 
critical  at  the  hump  and  sub-hump  speeds.  From  this  point  of  view,  the 
elbow  between  the  inlet  and  the  strut  is  the  critical  station  in  the 
system.  The  elbow  between  the  inlet  and  the  strut  necessarily  represents 
a radical  change  In  both  the  flow  direction  and  the  shape  of  the  passage. 
The  flow  is  turned  and  the  cross  section  of  the  duct  changes  from 
circular  to  a nearly  two-dimensional  form.  No  empirical  information  has 
been  found  on  the  performance  of  elbows  that  achieve  both  these  functions 

5 

simultaneously.  A study  of  available  data  on  turning  vane  corners 
recommended  thin,  noncircular,  arc  vanes  for  turning  the  flow.  A corner 
loss  coefficient  of  0.15  may  be  assumed  and  is  believed  to  be  attainable 
with  a reasonably  well  designed  elbow. 
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DESCRIPTION  OF  EXISTING  DESIGN  METHODS 


To  Illustrate  existing  inlet  design  capability,  this  section 
reviews  methods  which  have  been  used  to  generate  waterjet  inlet  designs. 
This  shall  elucidate  the  use  of  material  presented  in  the  preceding 
sections  and  will  hopefully  show  directions  in  which  improvements  can  be 
made. 

3 

Sherman  and  Lincoln  presented  a design  method  taking  into  account 
cavitation,  internal  and  external  losses.  The  authors  state  their  intent 
to  establish  a guideline  for  the  hydrodynamic  design  of  ram  inlets  and  to 
demonstrate  a logical  and  consistent  procedure  for  the  accounting  of 
external  drag  and  internal  system  losses.  Their  paper  does  offer  a very 
satisfying  explanation  of  design  philosophy,  based  on  sound  hydrodynamic 
principles.  The  occasional  unfortunate  abuses  of  syntax  gives  rise  to 
unnecessary  ambiguity. 

Certain  assumptions  are  made.  First,  it  is  assumed  that  an  optimum 
pump  flow  rate  has  been  established.  Second,  the  external  dimensions  of 
the  nacelle  and  strut  are  determined  by  cavitation  criteria:  the  absence 
or  presence  of  cavitation  for  the  subcavitating  and  supercavitating 
cases,  respectively.  For  a given  craft,  design  speed,  flowrate  and  IVR, 
the  external  nacelle  and  strut  dimensions  are  selected  according  to 
cavitation  requirements.  This  allows  little  latitude  in  designing  the 
internal  flow  path.  Therefore,  the  IVR  is  used  as  the  independent 
variable  to  arrive  at  an  optimum  external  and  internal  configuration.  By 
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Investigating  a range  of  IVR's  using  a cavitation  criterion,  they  are 
able  to  define  a particular  Inlet  system  for  each  case,  and  to  calculate 
Internal  and  external  performance  for  each  design,  thus  determining  the 
optimum  IVR  for  the  desired  speed  and  flowrate. 

The  selection  of  the  design  IVR  Is  shown  to  be  very  critical, 
making  It  necessary  to  adhere  to  the  optimum  value,  or  suffer  high  power 
losses.  At  100  knot  speeds,  base  vented  Inlets  appear  definitely 
superior  with  respect  to  drag  over  both  subcavltating  and  supercavitatlng 
types . 

The  Douglas- Neumann  program  Is  used  to  compute  the  potential  flow 
In  and  around  the  Inlet.  Considerable  leakage  1s  experienced,  a common 
observation  by  users  of  that  method.  From  these  computations,  the  mass 
flowrate  can  be  fixed  and  the  surface  pressure  distribution  determined. 

For  both  high  speed  subcavitating  and  base  vented  inlets,  the 
avoidance  of  cavitation  is  the  prime  design  criterion.  Cavitation  will 
affect  the  contour  of  the  nacelle  and  its  thickness  to  chord  ratio,  and 
the  shape  of  the  inlet  lip.  Most  important,  cavitation  will  influence 
the  selection  of  design  IVR  and  the  range  of  velocities  through  which  the 
inlet  is  able  to  operate  without  geometry  changes. 

As  a first  step  in  the  optimization  process,  values  for  the  inlet 
drag  coeffiqient  and  the  inlet  loss  coefficient  are  assumed.  A trade-off 
is  made  by  varying  the  flowrate  through  the  system  at  constant  power  and 
craft  speed.  The  optimum  system  is  selected  by  maximizing  net  thrust  or 
payload  capability  (reducing  system  weight).  These  calculations  are 
sunmarlzed  In  Figure  5. 
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et  Loss  Coefficient  = 0.2,  Design  IVR  = 0.75,  Design  Speed  = 50  knots 


1 


A computer  program  was  developed  which  would  take  as  Its  Input  the 
results  of  the  Neumann  program,  and  would  calculate  the  Inlet  performance 

3 

at  each  of  the  IVR's  of  interest.  This  program  gives  total  external 
drag  coefficients,  inlet  lip  cavitation  characteristics,  and  internal 
system  losses.  The  system  was  optimized  by  minimizing  the  combined  sum 
of  the  external  drag  losses  and  the  internal  flow  losses.  From 
investigations  with  the  Neumann  program,  it  was  ascertained  that  the  most 
satisfactory  inlet  shape  could  be  obtained  with  the  NACA  16-series 
forebody  thickness  distribution  generated  on  a camber  line  corresponding 
to  a surface  given  by 


(|  f * 

^0 


(^) 

^0 


2.4 


= 1 


where 


^0  0.793 


The  equation  for  the  forebody  of  the  NACA  16-series  is 


^ = 0.9897  (|)^^  - 0.2392  (|)  - 0.0410  (|)^  - 0.5594  (|)^ 
m 

Here,  t and  c can  be  used  to  scale  the  model, 
m 

The  required  thickness  to  chord  ratio  of  a NACA  16-series  section 
strut  was  calculated  from  the  following  formulas  for  a fully  wetted 
strut: 


t -1  + /I  + a. 
m _ 1^ 

‘^o  1.15 

where  is  the  incipient  cavitation  number, 


! 


1 

1 
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f 

and  for  a base  vented  strut: 

p 

m 1 

I I rWET^ 

I 

I The  procedure  for  obtaining  an  "optimum  inlet"  from  a given  set  of 

i design  conditions  is,  in  the  Sherman  and  Lincoln  method,  to  run  a matrix 

of  cases,  each  case  differing  from  the  others  in  geometrical  details. 

The  cavitation  characteristics  of  each  case  are  then  plotted.  From  these 
results,  the  proper  forebody  thickness  distribution  and  lip  thickness 
were  selected. 

Two  specific  ram  inlet  designs  were  initially  contemplated,  one  for 
a 50  knot  design  speed  and  another  for  a 100  knot  design  speed.  A range 
of  design  IVR's  from  0.25  to  0.875  was  considered.  It  was  decided  to 
concentrate  on  the  100  knot  design  because  it  was  felt  that  more  could  be 
learned  from  this  problem.  This  is  because  the  100  knot  design  requires 
more  area  for  external  diffusion  (to  keep  the  external  maximum  velocity 
within  given  limits),  a larger  variation  in  operational  IVR's,  and 
permitted  less  variation  in  local  velocity  above  freestream.  The 
operational  IVR's  for  the  two  designs  are  given  as  a function  of  craft 
speed  in  Figure  6. 

For  the  100  knot  design  speed  inlet,  an  optimum  flowrate  of  23,500  Ib/sec 
(10660  kg/sec)  was  determined.  Figure  7 shows  the  lip  contours  and 
Figure  8 shows  the  relationship  between  actual  nacelle  sizes.  In  the 
performance  calculations,  a 2 foot  (0.61  m)  strut  was  assumeo. 
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VARIATION  IN  OPERATIONAL  INLET  VELOCITY 
RATIO  WITH  VEHICLE  SPEED  AT  A CONSTANT  FLOW  RATE 

Design  Speed  = 100  kt 
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Figure  6B 
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The  shapes  of  the  cavitation  buckets  are  shown  in  Figure  9 as  a 
function  of  local  velocity  and  IVR.  This  figure  determines  the  operating 
conditions  under  which  the  inlet  is  free  of  cavitation. 

Figure  lOA  shows  the  additive,  system  and  pressure  drag 

coefficients,  components  of  the  external  drag,  as  a function  of  the 

design  IVR  for  a design  speed  of  100  knots.  Figures  lOB,  C,  and  D are 

the  corresponding  plots  for  the  base  vented  (at  maximum  thickness),  base 

vented  (forward  of  the  maximum  thickness)  and  supercavitating  ram  inlet 

systems,  respectively.  The  total  external  drag  coefficients  for  these 

inlets  are  compared  in  Figure  11,  where  one  can  see  that  the  ventilated 

design  is  superior  in  terms  of  external  drag  at  100  knots. 

The  inlet  loss  coefficient  K.  = — ^ — , where  H,  is  the  sum  of  all  the 

^ v//2g  *- 

internal  losses,  expressed  as  pressure  heads,  is  used  as  a measure  of 
internal  performance.  The  total  inlet  loss  coefficient  is  plotted  as  a 
funtion  of  the  design  IVR  in  Figure  12.  A breakdown  of  internal  losses 
by  component  is  shown  in  Figure  13.  The  variation  in  component  head 
losses  with  flowrate  at  an  IVR  of  0.75  are  shown.  The  nacelle  elbow 


losses  are  considerably  higher  than  those  of  the  other  components,  and 
all  losses  are  seen  to  vary  in  almost  direct  proportion  with  the  square 
of  the  flowrate.  From  Figures  10  and  11,  we  learn  that  the  external 
drag  decreases  with  design  IVR,  while  Figures  12  and  13  show  that  the 
internal  losses  Increase  with  increasing  IVR.  Remembering  that  design 


IVR's  correspond  to  different  inlet  shapes  (designs)  which  were  chosen 
according  to  their  cavitation  properties,  we  are  now  in  a position  to 
effectuate  a compromise  between  internal  and  external  losses  to  determine 
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Figure  9 - Cavitation  Characteristics  of  the  Six  Inlet  Designs 
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COMPARISON  OF  THE  TOTAL  EXTERNAL  DRAG  COEFFICIENTS 
OF  THE  INLET  SYSTEMS  INVESTIGATED  vs 
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INLET  VELOCITY  RATIO,  Vj/V^ 

Figure  12  - Total  Inlet  Loss  Coefficient  vs  Design  Inlet  Velocity  Ratio 
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the  best  Inlet.  Sherman  and  Lincoln  define  the  optimum  waterjet  Inlet  as 
the  one  which  constitutes  the  best  compromise  between  external  drag  and 
Internal  losses  while  also  meeting  the  desired  cavitation  requirements  of 
the  system. 

The  power  dissipated  by  Internal  losses  Is  given  by 

» 9 Q \ . >s  p V„3  (!i)  k,A, 

0 

while  the  power  associated  with  external  drag  Is 

“ ''o  ■ ’“>  Vi 

where  C_,  Is  the  total  external  drag  coefficient.  Since  the  term  JjpV  Is 
u 0 

common  to  both  factors,  we  can  compare  the  effects  of  Inlet  drag  and 
losses  by  plotting  CpA^  and  K^A^. (IVR)  as  a function  of  the  design  IVR. 

This  approach  has  been  employed  to  determine  the  optimum  configuration 
for  a ram  Inlet  with  a design  craft  speed  of  100  knots  and  a pump  flow 
rate  of  23,500  pounds  per  second.  Figures  14A  through  D are  plots  of  the 
available  trade-offs  and  the  resulting  optimum  design  IVR  for  the  basic 
Inlet  types  considered.  It  Is  obvious  from  these  figures  that  the 

subcavltating  and  supercavltating  Inlets  are  unsatisfactory,  and  that  ; 

base  venting  results  In  Improved  Inlet  performance.  A resume  of  the 
optimization  computations  Is  shown  In  Table  1. 

The  drag  performance  of  an  Inlet  at  100  knots  Is  extremely 
sensitive  to  the  choice  of  a design  IVR.  Consider  for  example  the  Inlet 
which  Is  base  vented  forward  of  the  point  of  maximum  thickness,  the  best 

performer  In  this  study.  Choosing  a design  IVR  which  Is  one  tenth  i 
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OPTIMUM  DESIGN  INLET  VELOCITY  RATIO  FOR  INLET 
BASE  VENTED  AT  THE  POINT  OF  MAXIMUM  THICKNESS 


INLET  VELOCITY  RATIO,  Vj /V 
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Figure  14D 
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Power  Loss 


Inlet  Type 

Optimum  Inlet 
Velocity  Ratio 

Inlet  Loss 
Coefficient 

Inlet  Drag 
Coefficient 

Horsepower 

Kilowatts  ' 

Subcavitating 

0,735 

0.268 

0.220 

10,750 

8,016 

Base  Vented 
at  Max  Thickness 

0.745 

0.285 

0.153 

9,520 

7,099  1 

Base  Vented  Forward 
of  Max  Thickness 

0.700 

0.215 

0.138 

7,825 

5,835 

Supercavitating 

0.720 

0.245 

0.285 

12,400 

9,247 

TABLE  1 • 

- RESULTS  OF  DESIGN 

INLET  VELOCITY 

OPTIMIZATION 

greater  than  optimum  results  in  a 23.4%  increase  in  power  loss.  A design 
IVR  which  is  ten  percent  smaller  than  optimum  increases  the  power  losses 
by  19%.  There  appears  to  be  a very  definite  optimum  design  inlet 
velocity  ratio  for  a 100  knot  ram  inlet  system  and  designing  even 
slightly  off  this  optimum  results  in  large  power  losses.  If  the  optimum 
base  vented  inlet  and  its  corresponding  design  IVR  are  selected,  very 
little  latitude  in  craft  speed  is  allowed  and  in  particular,  the  inlet 
will  choke  due  to  cavitation  with  a 20  knot  drop  in  craft  speed.  The 
important  conclusion  to  be  drawn  from  this  is  that  this  inlet  could  never 
accelerate  from  rest  to  its  design  speed  because  of  choking.  Hence,  a 
variable  geometry  inlet  is  an  absolute  necessity  in  this  case. 

Developmental  Sciences,  Inc.  (DSI)  has  reported  on  the  analysis  and 
design  of  pod  inlet. Their  reports  are  not  very  descriptive  and  do 
little  to  explain  their  inlet  design  method. 

From  the  given  thrust  requirements  of  the  propulsion  system  and 
preliminary  pump  characteristic,  the  (minimum)  pressure  recovery  factor 
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I Is  estimated,  using  Bernoulli's  law.  As  viscous  losses  are  not  Included 

t 

In  this  calculation,  the  figure  Is  likely  to  be  low.  Then  the  design 
cruise  IVR  Is  selected  and  cruise  and  hump  Inlet  areas  are  determined. 

• Then,  cavitation  choking  of  the  flow  In  the  96  degree  turn  Into  the  strut 

Is  examined.  For  this  consideration,  the  hump  (but  foll-bome)  operating 

I 

condition  Is  the  most  critical.  The  criterion  which  is  used  Is  that 
there  be  enough  diffusion  to  balance  any  cavitation  Inducing  pressure 
drop  across  the  turning  vanes. 

The  pod  and  plug  geometries  are  determined  primarily  by  the  duct 
area  and  the  Inlet  area  at  hump  and  cruise  conditions.  Where  the  plug  is 
fully  retracted,  the  Inlet  and  duct  exit  areas  are  required  to  be 
Identical.  It  Is  further  specified  that  the  open  area  distribution  along 

[ 

the  pod  between  these  two  points  be  constant.  This  implies  that,  (1) 
once  the  internal  area  distribution  of  the  pod  is  specified,  the  plug 
shape  is  determined,  (2)  the  shaped  portion  of  the  plug  must  have  the 
i same  length  as  the  diffuser  part  of  the  pod. 

I Therefore,  the  internal  shape  of  the  pod  surface  is  determined  by 

the  external  pod  shape  which  is  in  turn,  designed  on  the  basis  of 

i 

j cavitation-free  external  flow  at  100  knots.  The  diffuser  length  is 

I ideally  determined  by  the  requirement  that  the  flow  be  diffused  without 

I 

[ separation.  This  Is  considered  impractical  since  In  order  to  achieve 

fully  attached  flow,  the  pod  (and  hence  the  extended  centerbody)  become 
I unacceptably  long,  resulting  In  high  drag  and  structural  difficulties. 

Furthermore,  a longer  diffuser  means  a longer  plug  which  results  in 
greater  boundary  layer  development,  which  requires  an  even  longer 
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diffuser  to  avoid  separation.  Whether  these  statements  were  Investigated 
or  were  simply  speculations.  Is  not  mentioned. 


A balance  must  be  struck  between  short  and  long  diffusers  so  that 
at  the  cruise  deployment,  the  separation  point  Is  located  such  that  the 
required  total  recovery  Is  achieved. 

A number  of  Inlet  shapes  having  different  diffuser  lengths, 
contours,  etc.  were  generated  for  preliminary  screening  analysis.  By 
trial  and  error,  the  one  best  satisfying  the  requirements  was  selected 
and  analyzed  In  detail.  The  geometric  construction  of  the  final  selected 
Inlet  Is  shown  In  Figure  15.  A circular  arc  Is  used  to  describe  the 
external  pod  surface.  The  equation 

, V 144  r K 

<-nfc>  = * T7i?k-«7  > J 

where  k * 47.21637  and  = 20 

represents  the  Internal  pod  surface.  Ellipses  are  used  to  fair  the  inlet 
nose  Into  the  diffuser  and  external  pod  surfaces. 

A one-fifth  scale  model  of  the  DSI-desIgned  inlet  was  built  and 
tested  at  the  36  Inch  VPWT  at  DTNSRDC.®  The  pressure  recovery 
performance  of  the  inlet-diffuser  is  presented  as  a loss  coefficient  in 
Figure  16.  At  the  design  IVR  of  0.85,  the  pressure  loss  coefficient  is 
approximately  0.25,  which  compares  with  the  predicted  value  of  0.242  for 
the  cruise  condition.  For  the  low  speed  condition  corresponding  to 
simulated  speed  of  35  knots  and  a fully  retracted  centerbody,  the 
pressure  loss  coefficient  Is  shown  In  Figure  17.  Rather  high  losses  are 
Indicated.  At  the  hump  speed  design  IVR  of  .89,  the  value  of  the  model 
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PRESSURE  LOSS  COEFFICIENT 
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Figure  16-  Internal  Pressure  Loss  Performance  - High  Speed  Conditions 
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Figure  17  - Internal  Pressure  Loss  Performance  - Low  Speed  Conditions 


pressure  loss  coefficient  is  1.60.  The  test  report  attributes  this  to 
internal  cavitating  conditions.  This  value  greatly  exceeds  the  predicted 
value  of  0.282.  From  this  data,  it  is  concluded  that  the  inlet  design 
will  not  provide  the  required  performance  at  the  hump  speed. 

The  inlet  cavitation  inception  boundaries  for  the  cruise 
configuration  (centerbody  fully  extended)  are  shown  in  Figure  18. 

Cavitation  was  found  to  occur  on  the  inner  and  outer  surfaces  of  the 
inlet  lip  near  the  leading  edge.  No  cavitation  was  observed  on  the 
centerbody.  Therefore,  the  design  does  meet  its  requirements  at  100 
knots.  Furthermore,  the  experimental  data  and  prediction  for  external 
cavitation  are  in  good  agreement.  The  data  for  internal  cavitation 
indicate  a higher  attainable  IVR  than  predicted.  At  simulated  high 
speeds,  80  knots  or  better,  cavitation  was  observed  at  the  strut-pod 
intersection,  indicating  that  further  refinements  of  the  fairing  are 
required  in  this  region.  The  inlet  cavitation  inception  boundaries  for 
the  centerbody  fully  retracted  are  presented  in  Figure  19.  Internal  and 
external  cavitation  did  occur  on  the  inlet  lip  leading  edge.  At 
conditions  simulating  35  knots  - 18  feet  (5.49  m)  submersion,  cavitation  was 
observed  with  an  IVR  of  2.32  and  maximum  choked  flow  at  an  IVR  of  2.42 
(based  on  cruise  inlet  as  area).  At  simulated  35  knots  - 4 feet  (1.22  m), 
submergence,  internal  cavitation  inception  occurred  at  an  IVR  of  2.15  and 
choked  flow  at  an  IVR  of  2.20,  indicating  that  the  IVR  of  2.35  required 
to  accelerate  the  craft  through  the  hump  would  probably  not  be  achievable 
as  the  craft  rises. 
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Figure  18  - Inlet  Cavitation  Boundaries  Measured  and  Predicted, 
for  the  34.73  Inch  (0.882  tn)  Centerbody  Extension 
(Fully  Extended) 
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Figure  20  shows  typical  inlet  drag  performance  for  the  inlet  with 
fully  retracted  centerbody  at  a cavitation  number  of  2.83  simulating  20 
knots  full  scale.  Figure  21  shows  similar  data  for  the  fully  extended 
centerbody  at  a cavitation  number  of  0.082,  simulating  100  knots  full 
scale.  The  quantities  shown  are  the  total  measured  drag,  the  inlet 
momentum  drag, here  defined  as 
P q V. 

’ 

and  the  computed  external  drag  (»  measured  drag  - momentum  drag). 

5 ”9 

The  SES  design  manual  utilizes  the  method  of  Kuchemann  and  Weber 

to  determine  nacelle  dimensions.  A series  of  expression  is  developed  for 

bodies  with  different  configurations,  relating  the  maximum  external 


velocity  on  the  body  to  the  geometric  characteristics  of  the  nacelle  and 
inlet.  For  an  inlet  with  a fully  wetted  finite  afterbody,  the  formula  is 


DRAG  COEFFICIENT  - Cd 
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Figure  20  - Typical  Inlet  Drag  Performance  for  a Range  of  IVR, 
Centerbody  Retracted,  Simulated  20  Knots 
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Figure  21  - Typical  Inlet  Drag  Performance  for  a Range  of  IVR, 
Centerbody  Extended,  Simulated  100  Knots 


48 


where  o.  1s  the  incipient  cavitation  number,  A * the  maximum 
1 m 

cross-section  area,  A^.  Is  the  Inlet  area,  and  n Is  between  2 (for  long 

bodies)  and  3 (short  bodies).  Nothing  Is  said  about  specific  Inlet 

shapes,  and  drag  computations  are  not  mentioned. 

2 

The  basic  premise  In  the  Lockheed  Inlet  design  method  Is  that 
there  exists  an  optimum  relationship  between  Internal  system  losses, 
external  drag,  and  system  weight  such  that  the  propulsive  power  required 
to  produce  a given  thrust  Is  a minimum. 

The  Ideal  jet  propulsive  efficiency  1s^® 

_ work  per  unit  time  by  thrust 
^ ~ work  per  unit  time  by  ideal  pump’ 

Assuming  an  IVR  of  unity,  we  obtain 

p Q (Vj  - v„) 
n = “ 

’s  p Q (Vj^  - Vq^) 

where  p Is  the  fluid  density,  Q is  the  volume  flowrate,  Is  the  Inlet 
speed  (=  craft  speed  at  unit  IVR)  and  V.  Is  the  jet  speed.  This 
expression  reduces  to 


''“T-nTv; 

This  equation  Indicates  that  the  propulsive  efficiency  Is  Increased 
with  lower  values  of  jet  velocity,  which  corresponds  to  higher  values  of 
flowrates  for  a given  thrust.  However,  high  flowrates  require  a ducting 
system  with  large  cross-section  areas  to  minimize  the  flow  velocity. 

This  large  duct  result  In  prohibitive  values  of  both  external  drag  and 
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system  weight.  Therefore,  it  is  necessary  to  balance  the  internal  flow 
losses  of  a given  system  configuration  against  the  associated  external 
hydrodynamic  drag,  while  taking  into  account  the  resulting  propulsive 
efficiency.  Two  optimization  procedures  are  considered.  One  produces  a 
ship  with  the  minimum  shaft  horsepower  for  the  displacement,  speed  and 
drag  prescribed.  In  the  other,  the  system  is  optimized  for  the  maximum 
payload  with  a given  maximum  displacement. 

The  Lockheed  optimization  procedure  relies  heavily  on  drag 
estimates.  For  preliminary  design  purposes,  total  vehicle  drag  estimates 
are  obtained  first.  These  are  used  to  determine  power  requirements, 
engine  selection,  etc.,  and  to  obtain  rough  layouts  for  all  system 
components.  These  initial  drag  estimates  are  made  without  regard  to  the 
optimized  propulsion  system.  To  optimize  the  system,  subsequent  external 
drag  calculations  must  be  made  to  obtain  the  proper  trade-offs  between 
internal  and  external  flows.  The  resulting  shapes  and  drag  coefficients 
can  then  be  used  to  size  the  configuration,  consistent  with  the  system 
optimization  requirements.  Additional  external  drag  calculations  may  be 
necessary  to  optimize  a particular  component,  such  as  a diffuser,  again 
independent  of  the  system  optimization. 

The  relative  performance  of  selected  inlet  lip  profiles  can  be 
evaluated  on  the  basis  of  pressure  distributions  for  a range  of  inlet 
IVR's.  The  minimum  allowable  cavitation  number  is  determined  from  the 
operational  conditions  and  requirements  of  the  craft  and  the  internal 
flow  requirements  to  provide  the  necessary  thrust  over  the  entire  speed 
range.  Waterjet  inlets  should  be  designed  for  minimum  external  drag 
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within  the  limitations  of  internal  flow  and  cavitation  requirements. 

Minimum  drag  dictates  a small  external  diameter,  a minimum  lip  thickness, 
and  minimum  external  length.  Inlet  lips  are  relatively  short,  and  their 
length  is  normally  not  important  for  drag  considerations. 

A combination  of  internal  and  external  lip  shapes  has  been 
developed  numerically  to  satisfy  the  requirements  of  hydrofoil  vehicles. 

The  basic  precedure  used  is  the  numerical  solution  of  the  Neumann  problem 
for  potential  flow.  Once  the  basic  profiles  are  determined  by  trial  and 
error,  the  internal  and  external  shapes  are  geometrically  expanded  to 
provide  data  on  a family  of  inlets  with  varying  diameter,  length  and 
thickness.  Final  designs  for  actual  craft  may  require  further 
considerations  and  possibly  a final  computation  for  the  selected  inlet. 

The  Lockheed  design  method,  therefore,  relies  heavily  on  interaction  with 
the  designer. 

One  of  the  most  important  factors  to  be  considered  in  the  design  of 
high  speed  subcavitating  inlets  for  waterjet  propulsion  is  cavitation  of 
the  nacelle.  The  proper  IVR  must  be  selected  such  that  no  cavitation 
exists  for  the  desired  design  conditions  and  preferably  over  the  entire 
operating  range.  The  occurrence  of  cavitation  on  an  inlet  can  be 
controlled  by  the  proper  nacelle  contour  and  IVR.  For  increasingly 
greater  forward  speeds,  the  critical  pressure  at  which  cavitation  occurs 
decreased,  with  the  associated  requirement  for  more  streamlined  shapes. 

These  shapes,  however,  result  in  longer  nacelles  and  consequently  greater 
drag. 

j 

I 
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The  series  A Inlet  (see  Reference  9)  was  selected  as  the  basic 
profile  for  the  parametric  computer  calculations.  Several  members  of 
this  series  were  analyzed  with  the  Neumann  program  to  determine  the 
pressure  coefficient  distribution.  The  selection  of  a particular  design 
begins  with  the  determination  of  the  allowed  pressure  coefficient  from 
cavitation  requirements  and  the  known  operating  conditions.  From  the 
propulsion  requirements,  an  IVR  can  be  estimated.  With  this  IVR  and  the 
help  of  external  pressure  distribution  curves,  the  designer  selects  an 
inlet  size  which  will  avoid  cavitation,  i.e.,  he/she  chooses  a particular 
member  of  the  series  A on  the  basis  of  external  cavitation  performance. 
Once  this  shape  is  selected,  the  pressures  on  the  internal  lip  at  other 
operating  conditions  can  be  checked.  If  these  conditions  are  not 
satisfied,  it  may  be  necessary  to  select  another  shape  and  re-check. 

When  the  final  selection  is  made,  the  nacelle  size  is  determined.  The 
drag  must  then  be  calculated  and  checked  with  the  preliminary  design 
assumptions  to  determine  if  the  total  vehicle  drag  will  change 
significantly.  If  so,  another  iteration  may  be  necessary. 

Internal  and  external  flow  tests  were  made  on  a one  tenth  scale 
model  of  the  Lockheed  ram  inlet  design  for  a 550  ton,  80  knot 
supercavitating  hydrofoil,  and  for  a typical  50  knot  subcavitating 
hydrofoil.  The  models  represented  the  nacelle- inlet  foil -strut 
configuration  to  scale,  both  externally  and  internally,  up  to  the  pump. 

Subcavitating  model: 


The  internal  ducting  for  this  model  consists  of  a straight  wall 
conical  diffuser  having  a total  included  angle  of  12  degrees  and  an 


I 

I ^ (area)  diffusion  ratio  of  1.40  which  reduces  the  flow  velocity 

approximately  40%  prior  to  entering  the  elbow.  The  area  In  the  elbow  Is 

I 

I * then  reduced  by  20%  to  minimize  pressure  losses  and  match  the  available 

I 

area  in  the  strut.  The  strut  flow  is  then  diffused  with  a total  angle  of 
3.5  degrees  and  an  area  ratio  of  1.90.  The  strut  flow  Is  divided  into 
^ three  channels  by  two  splitters,  the  location  of  which  is  such  that  the 

I diffusion  rate  is  the  same  for  each  channel.  The  splitters  also  help 

provide  structural  rigidity. 

Towing  basin  tests  were  run.  At  the  scaled  cruise  velocity  of  50 
knots,  the  model  was  found  to  be  free  of  cavitation  at  IVR's  of  .52  and 
.70.  Similarly,  no  cavitation  was  found  at  30  knots,  the  take-off  speed. 

The  overall  pressure  loss  coefficient,  excluding  elevation  losses,  is 
shown  Is  Figure  22,  plotted  against  the  Reynolds  number  based  on  inlet 
diameter  and  velocity.  Total  inlet  drag  is  presented  in  Figure  23  as  a 

1 

function  of  IVR  for  two  values  of  the  submergence.  Interestingly,  the 
foil  accounts  for  about  40%  of  the  total  drag. 

1 
1 

Supercavitating  model: 

Provisions  were  made  to  allow  the  testing  of  various  configurations 
of  the  system  by  interchanging  the  following  duct  components: 

(a)  two  inlet  lip  contours, 

(b)  three  elbows  with  varying  rake  angles, 

(c)  a variety  of  turning  vanes,  and 

(d)  two  strut  diffusion  ratios. 
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DRAG  COEFFICIENT 


This  model  was  equipped  with  supercavitating  foils.  The  nacelle  itself 
was  not  meant  to  be  operated  under  cavitating  conditions.  At  76  knots 
with  an  IVR  of  0.52  and  a submergence  of  8.4  feet  (2.56  m),  the  nacelle  was 
observed  to  cavitate.  In  general,  cavitation  could  be  delayed  by  either 
Increasing  the  IVR  or  the  submergence.  Thus  increasing  the  IVR  to  0.78 
eliminated  cavitation  In  this  case.  The  cavitation  margin  for  these 
conditions  was  narrow.  Indicating  that  further  speed  Increases  with  this 
inlet  are  likely  to  give  rise  to  cavitation.  It  is  doubtful  whether  this 
inlet  could  be  successfully  operated  at  100  knots. 

The  overall  pressure  loss  coefficient  for  the  supercavitating  model 
includes  all  losses,  with  the  exception  of  cavitation,  from  free  stream 
to  strut  exit,  and  are  plotted  in  Figure  24  against  Reynolds  number.  In 
addition,  the  loss  coefficient  to  the  pump  is  also  indicated.  The  drag 
coefficients  obtained  for  the  supercavitating  model  are  shown  in  Figure 
25. 

RAM  INLET  PERFORMANCE 

The  subject  of  ram  inlet  performance  is  a particularly  difficult 
one  to  address  at  this  time  in  view  of  the  (surprising)  paucity  of 
experimental  data.  The  performance  data  presented  in  this  section  are 
derived  exclusively  from  experimental  results.  This,  it  was  felt,  could 
provide  the  most  reliable  way  of  comparing  various  concepts  for  the 
purpose  of  singling  out  a particular  concept  on  which  to  concentrate 
further  research  and  developmental  efforts.  Three  inlets  are  compared  in 
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A OVERALL  LOSS  COEFFICIENT  TO  CARRIAGE  PUMP  INLET 
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Figure  25  - Drag  Coefficient  for  Supercavi fating 
Lockheed  Model 
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this  section:  the  DSI  pod  with  moveable  centerbody,  the  Lockheed  inlet, 
and  the  Aerojet  inlet  used  on  the  SES  lOOA.  The  data  for  these  inlets 
were  taken  primarily  from  References  2,  7,  8,  11,  and  12. 

The  successful  operation  of  strut-pod-propelled  craft  such  as  the 
Aerojet  SES  lOOA  and  the  PGH-2  Tucumcari  has  demonstrated  the 
practicality  of  such  inlets.  The  design  technology  for  ram  inlets  is, 
in  view  of  all  the  available  experimental  evidence 
obtained  thus  far,  aooroximatelv  as  advanced  as  that  for 
flush  inlets.  The  Aerojet  SES  lOOA  has  achieved  a maximum 
speed  approaching  the  design  maximum  speed  with  waterjet  propulsion  and  a 
base-vented  strut-pod  inlet.  The  strut  and  pod  are  base-ventilated  to 
achieve  cavitation-free  performance  at  the  high  maximum  speed  and  to 
reduce  inlet  system  drag.  An  auxiliary  inlet  located  in  the  strut 
leading  edge  provides  variable  inlet  area. 

The  SES  lOOA  model  test  data,  while  limited,  are  useful  in  setting 
an  upper  bound  on  inlet  drag  for  very  high  speed,  variable  area  inlets. 
These  data  indicate  a high  speed  inlet  drag  coefficient  (based  on  high 
speed  inlet  area  and  ship  speed)  of  about  0.45.  The  energy  recovery  data 
are  complicated  by  the  complex  internal  ducting.  The  data  of  Reference 
13,  for  a simplified  geometry,  indicate  a loss  coefficient  of  0.10  to 
0.15  at  a typical  design  IVR  of  0.7-0. 8 for  the  initial  portion  of  the 
inlet.  From  these  data,  overall  loss  coefficients  of  0.20  to  0.25  seem 
probable  at  typical  dfesign  IVR's.^^ 

Three  inlet  lips  were  tested  with  inlet-to-pod  diameter  ratios  of 
0.4,  0.45  and  0.50.  For  each  of  these,  the  internal  and  external 
cavitation  boundaries  were  determined.  Figure  26  indicates  the 
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cavitation  performance  of  each  diameter  ratio  inlet.  Based  on  its  high 

cavitation- free  internal  flow  capability,  the  0.50  diameter  ratio  inlet 

was  selected  for  complete  performance  evaluation.  The  auxiliary  inlet 

was  tested  at  the  full  open  (8.5  inches,  .216  m)  position  in  conjunction  with  the 

0.45  diameter  ratio  primary  inlet.  Within  the  speed  and  IVR  ranges 

investigated,  the  primary  inlet  would  cavitate  internally  before 

inception  took  place  on  the  auxiliary  inlet,  see  Figure  27. 

The  cavitation  performance  of  the  DSI  pod  was  found  to  be 
satisfactory  at  scaled  cruise  speed,  but  the  inlet  choked  at  scaled  hump 
speed.  From  this,  it  was  concluded  that  the  inlet  could  not  accelerate 
the  craft  beyond  the  hump  condition.  Although  it  has  been  speculated 
that  a redesign  of  the  centerbody  could  solve  this  problem,  no  steps, 
theoretical  or  experimental,  have  been  taken  in  this  direction. 

The  drag  figures  reported  for  various  inlet  concepts  are  not  easy 
to  interpret.  For  example,  the  drag  coefficients  which  are  presented  for 

q 

the  DSI  pod  are  based  on  estimated  total  wetted  surface  area.  Drag 
figures  obtained  from  different  tests  cannot  be  compared  since  typically, 

q 

different  inlet  components  are  included.  For  example,  Sobolewski 

reports  on  measurements  performed  on  the  pod  plus  a nearly  vertical 

13 

strut,  while  the  test  performed  on  the  Aerojet  inlet  included  the  pump 
and  nozzle.  While  shaft  and  nozzle  thrust  were  "measured  out,"  pump, 
ducting  and  exit  nozzle  drag  were  not  included  in  the  final  drag  figures. 

The  effects  of  interactions  near  the  impeller  can  only  be  guessed  at. 

For  these  reasons,  it  is  not  possible  to  meaningfully  compare  drags  from 
various  tested  inlet  concepts.  Similarly,  the  lack  of  uniformity  in 
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CAVITATION  SCALED  SHIP  SPEED  - KNOTS 


reporting  pressure  losses  makes  It  difficult  to  establish  comparisons. 
Cavitation,  on  the  other  hand,  Is  less  ambiguously  reported,  simply 
because  of  the  qualitative  nature  of  the  data.  A performance  chart.  Table 
2,  comparing  the  three  tested  Inlet  systems  summarized  the  relevant 
experimental  data.  On  the  basis  of  cavitation  alone,  the  Aerojet  inlet  Is 
the  only  one  which  performed  satisfactorily;  presumably,  the  DSI  and 
Lockheed  Inlets  could  be  redesigned  to  avoid  cavitation,  but  the  Aerojet 
Inlet  Is  already  full-scale  proven  concept. 

The  most  Important  performance  criterion  for  waterjet  Inlets  Is  that 
they  operate  without  choking  at  hump  condition  and  without  extensive 
cavitation  at  cruise  speed.  An  Inlet's  first  requirement  Is  that  it 
operates  without  unwanted  cavitation.  Secondary  requirements  Include 
efficient  performance  (low  drag,  good  pressure  recovery).  From  this  point 
of  view  the  Aerojet  Inlet  Is  superior  to  Its  candidates,  with  a 
respectable  cavitation  margin  aIVR  = 0.2  at  100  kt  with  an  IVR  of  0.70. 

The  Lockheed  Inlet  system,  which  barely  avoids  cavitation  at  76  kt  with  an 
IVR  of  0.75,  and  the  DSI  Inlet  with  Its  hump  cavitation  choking,  are 
poorer  performers.  It  would  appear  that  Aerojet's  strut  auxiliary 
openings  constitute  an  effective  and  proven  method  of  Increasing  the 
flowrate  at  hump  condition. 
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DISCUSSION  AND  RECOMMENDATIONS  FOR  INLET  DESIGN 
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Three  ram  Inlet  design  methods  were  reviewed.  The  DSI  method 
designs  an  inlet  according  to  cavitation  requirements  and  then  calculates 
the  resulting  drag.  Should  the  inlet  be  too  draggy,  modifications  are 
made  on  the  basic  design.  The  Lockheed  method  is  similarly  not  suited  for 
computer  implementation  as  it  relies  heavily  on  humans  flipping  through 
graphs  and  charts.  The  Sherman  and  Lincoln  method,  however,  can  be 
adapted  for  numerical  optimization.  Here,  the  various  drag  formulas  can 
be  entered  as  separate  program  elements  and  the  optimization  can 
effectively  be  made  by  computing  the  trade-offs  between  internal  losses 
and  external  drag. 

Certain  incertitudes  remain,  especially  with  regard  to  the  relation 
between  pressure  and  additive  drag.  Clarifications  should  be  made  in  this 
area.  In  addition,  the  drag  computation  method  could  be  updated  to 
include  recent  developments  in  this  field.  With  these  improvements,  the 
Sherman  and  Lincoln  method  could  be  adapted  as  the  Navy's  in-house  design 
method  for  ram  inlets. 

14 

According  to  Etter,  for  speeds  less  than  50  knots,  fixed  geometry 
waterjet  inlets  are  acceptable.  In  the  case  of  high  performance  craft, 
especially  those  with  low  speed  humps,  the  inlet  design  may  be  governed  by 
the  need  for  variable  geometry  in  order  to  accommodate  a cavitation-free 
passage  through  the  hump  speed.  For  SES  and  hydrofoil  ships,  which  have 
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nearly  equal  flow  requirements  at  hump  and  cruise  speeds,  the  inlet  must 

either  have  a variable  area  capability  or  be  able  to  operate  over  a wide 

range  of  IVR's.  Etter  maintains  that  the  most  practical  variable  area 
scheme  for  pod-strut  inlets  is  the  use  of  an  auxiliary  strut  inlet  which 

can  be  opened  or  closed  with  a sliding  gate.  The  only  other  variable  area 

concept  which  appears  feasible  for  high  speed  pod-strut  inlets  is  the 
"sliding  centerbody"  concept. 

Barr  and  Stark^^  reviewed  the  performance  of  existing 
water jet-propelled  high  speed  craft.  The  SES  lOOA  utilizes  base-vented 
strut-pod  inlet.  Variable  geometry  is  achieved  by  means  of  an  auxiliary 
inlet  located  in  the  strut  leading  edge,  allowing  a 135%  area  increase 
over  that  of  the  pod  alone. 

For  high  speed  craft  such  as  hydrofoils  and  SES,  which  have  large 
thrust  requirements  at  low  speed,  the  use  of  variable  area  inlets  is 
essential.  As  the  flowrate  requirements  at  hump  and  cruise  speeds  are 
approximately  equal  in  the  applications,  a wide  variation  in  operating 
IVR's  is  needed  to  accommodate  the  flow.  A fixed  area  inlet  operating  over 
a wide  range  of  IVR's  would  necessarily  have  a large  inlet  radius,  which 
in  turn  significantly  adds  to  the  inlet  drag.  A completely  variable 
inlet,  on  the  other  hand,  can  accommodate  a constant  flowrate  over  the 
entire  speed  range  without  any  change  in  IVR,  and  requires  no  leading  edge 
radius.  Such  an  inlet  has  a lower  drag  coefficient.  A realistic  design 
is  a compromise  between  these  two  extremes,  with  a finite  leading  edge 

radius  allowir\g  operation  over  a limited  range  of  IVR's. 

14 

In  a review  of  waterjet  propulsion,  Etter  remarks  that  the  strut 
cross  section  must  provide  adequate  structural  strength  and  internal  duct 


area  and  must  permit  cavitation- free  operation  at  all  operating  speeds  and 
yaw  angles.  Below  60  knots,  fully  wetted  sections  have  the  lowest  drag. 

At  high  speeds,  fully  wetted  sections  have  poor  strength  and  Inadequate 
resistance  to  cavitation  at  yawed  conditions  due  to  their  small  thickness, 
making  base-vented  sections  preferable. 

The  comparison  of  various  design  concepts  Is  a difficult  task  for 

several  reasons.  Inlet  design  Is  still  In  Its  infancy  and  the  current 

state-of-the-art  leaves  plenty  of  room  for  improvements,  particularly  in 

the  area  of  drag  and  cavitation  predictions  with  flow  over  complex 

boundaries.  The  comparison  of  experimental  data  is  not  so  clearly  defined 

either.  For  example,  Boeing  designed  an  Inlet  system  for  the  PGH-2. 

These  fixed  area,  rectangular-shaped  Inlets  were  designed  for  a cruise 

speed  of  50  knots  and  tested  at  the  St.  Anthony  Falls  Hydraulic 
15 

Laboratory.  The  pressure  loss  coefficient  was  measured  over  IVR's 
ranging  from  1.20  to  2.5.  The  lack  of  data  below  unit  IVR  makes 
comparisons  difficult  with  other  data.  Furthermore,  the  P6H-2  inlet  has 
several  strut  elbows,  a lossy  duct  geometry.  On  the  other  hand,  the  craft 
operated  satisfactorily  for  several  years.  No  serious  cavitation  damage 
occurred. 

It  is  recommended  that  the  Sherman-Lincoln  design  philosophy  be 
utilized  In  developing  a computer  code  for  the  generation  of  axisymmetric 
strut  pod  Inlet  systems.  The  program  should  have  a modular  struture  to 
facilitate  Improvements (see  for  example  Reference  16).  The  nacelle  shape  families 
and  the  drag  estimation  methods  should  be  reviewed  in  the  light  of  recent 
developments. 

Based  on  admitedly  limited  test  data  on  auxiliary  inlets  and  variable 
geometry,  it  appears  that  the  high  speed  aerojet  concept  is  superior. 
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It  is  therefore  recommended  that  the  high-speed  Aerojet  concept  be 


adopted  for  further  development.  This  system  utilizes  an  auxiliary 
opening  situated  on  the  strut  leading  edge  and  operated  by  a sliding  gate. 
Unlike  the  movable  centerbody  concept,  no  increase  in  nacelle  size  is 
needed  to  accommodate  the  increased  flowrate  mechanism.  This  concept  was 
successful  in  the  SES  lOOA  and  therefore  has  proven  potential  as  a parent 
form  for  similar  strut  pod  inlet  systems.  The  preliminary  developmental 
design  work  performed  by  Code  1532  can  be 
applied  to  the  design  of  such  inlets. 

When  high  speed  (=  100  kt)  applications  are  considered,  Sherman  and 

3 

Lincoln  have  convincingly  demonstrated  the  need  for  a variable  inlet  for 

cavitation- free  operation.  For  relatively  low  speed  applications  (<  50 

kt)  where  a significant  hump  drag  does  not  exist,  fixed  area  waterjet 

14 

inlets  may  be  acceptable,  according  to  Etter.  In  any  case,  caution 
should  be  exercised  when  designing  inlets  for  the  upper  end  of  this  speed 
range,  where  additional  inlet  area  may  prove  beneficial  at  off-design 
conditions.  This  makes  the  Aerojet  auxiliary  strut  apertures  even  more 
attractive  since  they  represent  an  easier  design  problem  than  the  fitting 
of  a moveable  centerbody.  In  designing  for  this  speed  range,  a fixed  area 
inlet  should  first  be  considered.  If  adequate  flowrate  is  not  provided  at 
off-design  conditions,  then  auxiliary  inlets  can  be  added  to  the  design 
without  altering  the  nacelle  Itself.  Thus,  for  the  intermediate  speed 
range,  it  is  recommended  that  the  benefits  of  having  additional  intake 
area  be  investigated  on  a case-by-case  basis. 
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ANALYTICAL  DRAG  PREDICTIONS 


The  success  of  any  1n1et  design  method  to  a large  extent  depends  on 
Its  ability  to  estimate  drag  accurately.  This  section  reviews  drag 
calculations  within  the  context  of  waterjet  Inlet  design.  We  first 
present  a brief  discussion  of  the  various  forms  of  drag  (e.g.,  skin 
friction  drag,  form  drag,  etc.),  emphasizing  fundamental  hydrodynamic 
principles.  A brief  review  of  current  theoretical  predictive  techniques 
follows.  Design  programs  usually  utilize  drag  prediction  methods  which 
are  based  on  empirical  knowledge,  an  attempt  to  fill  in,  by  careful  usage 
of  experimental  data,  the  large  gaps  In  the  formal  theories.  Various 
design  programs  will  be  reviewed  and  their  drag  prediction  methods 
outlined. 

Skin  Friction  Drag 

Viscosity,  through  the  action  of  shear  force,  couples  the  body 
surface  to  the  fluid,  allowing  the  body's  momentum  to  diffuse  throughout, 
resulting  In  fluid  acceleration.  The  body  therefore  experiences  a force, 
called  skin  friction  or  viscous  drag,  contrary  to  the  direction  of  flow  In 
the  vicinity  of  Its  surface.  Where  the  flow  Is  laminar,  viscous  diffusion 
of  momentum  Is  the  only  source  of  skin  friction  drag,  and  a reliable 
theory  has  been  developed.  With  turbulent  flow,  eddies  contribute  to 
momentum  transfer  by  a convective  process  for  which  no  rigorous  theory 
exists  at  present.  However,  many  useable  empirical  results  exist  (e.g. 
the  Schoenherr  line). 
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Form  Drag 

Whereas  skin  friction  drag  is  (locally)  a tangential  force,  a normal 
force  of  viscous  origin  can  also  arise  when  flow  occurs  past  a surface 
which  is  not  everywhere  parallel  to  the  main  stream.  This  force  is  known 
as  the  pressure  drag  or,  since  it  depends  on  the  shape  of  the  boundary, 
the  "form"  drag.  Another  name  often  used  is  "viscous  pressure  drag," 
which  necessarily  excludes  the  induced  drag  due  to  lift. 

An  Important  cause  of  form  drag  is  flow  separation.  A boundary 
layer  in  an  adverse  (positive)  pressure  gradient  will  experience  a 
retarding  effect  which  will  be  particularly  strongly  felt  near  solid 
boundaries,  where  the  fluid  possesseslittle  momentum.  In  fact,  the  layer 
of  fluid  immediately  adjacent  to  the  solid  boundary  may  lose  all  of  its 
streanwise  momentum,  and  reverse  its  direction  under  the  influence  of  the 
adverse  pressure  gradient  . The  recirculation  (separation)  regions  thus 
created  are  typically  at  lower  pressures  than  we  would  expect  from  purely 
potential  flow  pressure  recovery.  Thus,  pressure  differences  can  be 
induced  across  a body  by  flow  separation,  resulting  in  a net  drag. 

At  the  base  of  projectiles,  there  originates  a pressure  drag  which 
is  consequently  termed  "base  drag."  This  is  a result  of  the  boundary 
layer  separation  at  the  base  edge.  A good  discussion  of  this  phenomenon 
is  given  by  Hoerner.^^ 

Cavity  Drag 

The  drag  of  bodies  that  are  exposed  to  cavitation  or  ventilation,  is 
not  the  same  as  that  for  the  fully  wetted  case.  Generally,  cavitation  is 


70 


1 


a condition  to  be  avoided  in  the  design  of  strut-pod  inlets.  There  are 
situations  however  when  ventilation  cannot  be  avoided,  or  when  cavitation 
(or  ventilation)  are  deliberately  introduced  to  reduce  the  wetted  area 
even  at  the  expense  of  increased  cavitation  drag.  It  is  with  the  last 
situations  in  mind  that  we  now  consider  the  subject  of  drag  at  low 
cavitation  numbers. 

The  cavity  drag  of  a body  corresponds  to  an  average  positive 

] 

pressure  on  its  wetted  side  and  to  the  uniform  negative  pressure  within 

the  cavity. The  drag  coefficient  of  sharp-edged  bodies  (e.g.,  discs  and 

cones)  can  be  reliably  predicted  with  the  relation 

Cp  = Cp  (1  + a) 

0 

where  a = (p^  - Pyapor^/^*  cavitation  index,  and  Cp  is  the  drag 

coefficient  for  the  a = 0 condition.  This  formula  is  also  valid  for 

IR 

supercavitating  wedge-shaped  struts. 

The  precise  effect  of  cavitation  number  on  the  drag  of  even  simple 
bodies  of  revolution  with  fixed  detachment  is  not  entirely  understood. 

The  drag  coefficient  is  known  to  increase  asymptotically  (a  -*■  0)  according 
to  the  above-given  formula.  However,  the  magnitude  of  the  a range  over 
which  this  law  is  applicable  depends  on  the  fullness  of  the  body. 

Interference  Drag 

This  is  an  inherently  "engineering"  term,  a powerful  reminder  of  the 
complexities  and  apparent  vagaries  of  fluid  dynamics,  and  a testimony  to 
the  engineer's  determination.  Interference  drag,  unlike  skin  friction. 
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[ form  drag  or  cavity  drag,  is  not  related  to  any  specific  physical 

I phenomenon.  Rather,  it  arises  when  complex  flows  are  regarded  as  being 

formed  of  a superposition  of  simple  (and  calculable)  flows,  plus  an  ad  hoc 
"interference"  or  interaction  term  to  account  for  any  difference  between 
this  simplistic  and  asynergistic  model  and  reality. 

A good  review  of  hydrodynamic  interactions  on  high  performance  craft 
propulsion  systems  is  given  in  Reference  19. 

Wave  drag,  the  resistance  associated  with  energy  being  radiated  away 
in  the  form  of  surface  gravity  waves,  is  usually  considered  negligible 
for  strut-pod  inlets,  since  wave  drag  decreases  markedly  with  increasing 
submergence.  Furthermore,  at  the  high  Froude  numbers  encountered  wave 
drag  is  expected  to  be  small, 

PREDICTIVE  STATE-OP-THE-ART 

In  this  section,  we  describe  analytical  methods  representative  of 
our  current  state  of  knowledge  in  the  area  of  hydrodynamic  drag 

prediction.  Rather  than  being  fully  understood,  the  major  contributing  j 

i 

factors  to  drag  are  subjects  of  active  research.  Few  are  the  instance  ! 

where  we  stand  on  firm  theoretical  ground  and  consequently,  empiricism  is 
often  used. 

Excluding  cavity  and  lift  effects,  hydrodynamic  drag  is  a 
consequence  of  fluid  viscosity.  The  relationship  between  viscosity  and 
drag  may  be  a direct  one,  as  in  the  case  of  skin  friction  drag,  or 
indirect,  as  where  flow  separation  or  turbulence  occur.  The  concept  of 
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the  boundary  layer,  either  laminar  or  turbulent,  is  almost  always  used  so 

that  potential  flow  calculation  methods  do  find  their  place  in 

hydrodynamic  drag  predictions,  as  a first  step  in  the  evaluation  of 

boundary  layer  effects.  Potential  flow  has  been  extensively  studied  and  a 

mathematically  rigorous  theory  has  been  firmly  established. 

Unfortunately,  the  small  number  of  practical  cases  for  which  analytic 

solutions  can  be  obtained  precludes  a general  analytical  approach  to 

design  problems.  With  the  advent  of  digital  computers  and  efficient 

numerical  techniques,  potential  flow  problems  with  complicated  boundaries 

became  tractable.  Two  general  approaches  can  be  taken.  In  one,  the 

equations  of  Laplace  (or  Poisson,  where  sources  are  present)  are  directly 

solved  by  finite  difference  schemes  within  boundaries  on  which  the  flow  is 

prescribed.  The  solution  must  be  found  at  every  mesh  point  even  though 

one  may  only  be  interested  in  the  flow  over  a certain  limited  subdomain, 

over  a body  surface  for  example.  This  method  will  not  be  considered 

here.  An  alternative  and  mure  economical  approach  consists  in  replacing 

bodies  with  appropriate  source  distributions  giving  the  same  normal  flow 

at  the  boundaries.  The  total  number  o'  nodal  points  is  therefore  reduced 

to  that  necessary  to  describe  body  surfaces,  and  the  problem  reduces  to  a 

computation  of  the  source  strength,  from  which  the  flow  at  any  point 

external  to  the  body  may  also  be  found.  This  approach  has  been  refined  by 

20 

J.C.  Hess  and  A.M.O.  Smith.  Axial  symmetry  can  be  exploited  to  further 
reduce  the  number  of  surface  elements  and  allowances  can  be  made  to 
introduce  in  the  computation  those  flow  components  which  can  be  described 
analytically,  resulting  in  an  increased  computational  flexibility.  Thus, 
potential  flow  calculations  are  well  in  hand. 
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A distinction  must  be  made  between  theoretical  calculations  used  as 
research  tools  and  design  computational  methods.  Rather  than  offering  an 
exhaustive  review  of  recent  advances  In  the  fundamental  research  of 
boundary  layers,  flow  stability  and  other  topics  of  significant  relevance 
to  drag  prediction,  we  choose  rather  to  concentrate  on  calculatlonal 
methods  actually  used  by  designers,  and  the  drag  models  used  to  motivate 
these  techniques;  we  pay  only  cursory  lip  service  to  research 
calculations. 

Reliable  methods  for  the  numerical  computation  of  boundary  layers 
have  been  developed. These  solve  the  boundary  layer  equations  1n 
their  differential  form  and  therefore  provide  information  about  details  of 
the  flow.  These  methods  can,  and  indeed  have  been  used  in  drag 
calculations,  but  require  numerous  computational  steps  resulting  in  much 
information,  little  of  which  is  directly  relevant  to  viscous  drag.  The 
reason  is  that  drag  is  an  integral  property  of  the  flow  and  can  be 
computed  without  a detailed  knowledge  of  the  boundary  layer  structure. 

More  economical  approaches  in  drag  predictions  have  evolved  from  so  called 
integral  methods  in  which  the  boundary  layer  equations  are  first 
integrated  in  the  cross-stream  direction  to  provide  relationships  between 
various  gross  features  of  the  flow. 

The  viscous  friction  experienced  by  a body  in  motion  in  a real  fluid 
can  be  obtained  from  the  body's  influence  on  the  surrounding  fluid. 

Through  the  action  of  viscosity,  momentum  is  imparted  to  the  fluid, 
evident  at  first  as  large  velocity  gradient  in  the  boundary  layer.  The 
effect  is  cummulative  as  one  proceeds  from  the  leading  edge  stagnation 
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point  to  downstream  stations.  Increments  in  body  surface  result  in 
increased  viscous  drag  and  additional  shed  momentum,  which  diffuses  and 
contributes  to  the  thickening  of  the  boundary  layer.  Clearly  then,  an 
evaluation  of  the  accummulated  momentum  deficit  at  the  trailing  edge  or  in 
the  wake  can  be  utilized  in  calculating  the  viscous  drag.  Such  methods 
have  been  developed  by  integral  calculations  and  have  produced  various 
"approximate  formulae"  relating  drag  to  body  shape,  potential  flow  and 
kinematic  viscosity.  As  a rule,  models  ignore  separation  and  are 
regime-dependent.  That  is,  a uniformly  valid  formula  for  laminar  and 
turbulent  flow  cannot  be  found.  Although  these  formulas  are  relatively 
accurate,  their  applicability  to  practical  cases  depends  upon  reliable 
methods  to  determine  the  point  separating  the  laminar  and  turbulent  flow 
regions.  Often,  laminar  separation  trigger  turbulence  so  that  its  onset 
must  also  be  considered. 

In  the  drag  prediction  methods  whose  descriptions  follow,  a pattern 
is  seen  to  emerge.  The  flow  is  divided  into  various  regions:  potential 
away  from  the  body,  a laminar  boundary  layer  grows  from  the  stagnation 
point  up  to  a transition  point  whose  location  is  determined  by  one  of 
several  semi-empirical  techniques.  Thereafter,  the  flow  is  assumed  to  be 
fully  turbulent. 

Laminar  separation  presents  additional  problems  which  have  been 
conveniently  dismissed  with  the  assumption  that  the  elimination  of 
separation  form  drag  represents  the  first  step  in  finding  optimal  body 
shapes.  Consequently,  most  design  methods  are  concerned  with  reducing  drag 
on  nonseparating  bodies. 
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To  be  sure,  laminar  separation  can  be  predicted  numerically,  and 
shows  up  as  a singularity  beyond  which  boundary  layer  computations  cannot 


be  continued.  The  singularity,  however  is  a consequence  of  the 
mathematical  formulation  used,  and  not  a feature  of  the  flow  to  be  studied. 
Nevertheless,  it  does  serve  as  an  indicator  of  separation.  The  parabolic 
nature  of  the  boundary  layer  equation  precludes  a continuation  of  the 
solution  within  the  separated  region  in  the  absence  of  upstream  boundary 
conditions.  These  and  other  topics  related  to  our  current  state  of 


knowledge  on  boundary  layer  are  lucidly  treated  in  two  excellent  review 
25,  26 


.27 


papers. 

J.L.  Hess  presented  a method*"'  to  design  by  analytic  means  a class  of 
axi symmetri c bodies  having  low  drag  in  incompressible  flow  for  the  case 
where  the  boundary  layer  is  fully  turbulent  over  the  entire  body.  He 


compared  various  drag  calculation  methods  to  find  that  the  Truckenbrodt 
28 


formula,  which  expresses  drag  as  an  integral  of  a power  of  the  potential 
surface  velocity,  is  sufficiently  accurate,  and  this  formula  is  adopted  as 
his  chief  analytical  tool: 

^ 10/3 
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where  the  Integral  is  over  the  entire  body  profile.  L is  the  body  length, 
u^  the  potential  surface  velocity,  s the  arc  length  and  is  the  skin 
friction  coefficient,  which  for  smooth  bodies  is  taken  to  be 


0.455 


(log  N^) 
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where  Np  1s  the  Reynolds  number  based  on  body  length.  This  expression 
differs  little  from  the  friction  coefficient  obtained  from  the  Schoenherr 
line  as  expressed  by  the  formula 


log  (Np  C^)  * 


0.242 


A comparison  of  the  two  appears  on  Table  3. 


log  Rg 

Schoenherr 

X 10^ 

455 

(log  Np)^'^® 

5 

7.18 

7.16 

6 

4.40 

4.47 

7 

2.93 

3.00 

8 

2.07 

2.13 

9 

1.53 

1 .57 

10 

1.17 

1.20 

TABLE  3 - SKIN  FRICTION  COEFFICIENTS 

In  this  paper,  three  drag  prediction  methods  were  compared  with 
experimental  data  to  determine  the  relative  accuracy  of  these  methods.  The 
Truckenbrodt  formula  was  compared  with  the  result  of  using  the  finite 
difference  boundary  layer  method  in  Young's  formula.  The  outcome  of  this 
comparison  was  that  the  simplified  integral  formula  of  Truckenbrodt  is 
definitely  more  accurate.  Next,  the  momentum  integral  boundary  layer 
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method^ Patel  ,was  used  in  conjunction  with  Granville's  drag  formula 
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This  was  compared  with  Truckenbrodt's  formula,  and  with  the  finite 
difference-augmented  Young  formula.  The  calculation  based  on  the  momentum 
integral  (Patel -Granville)  was  found  to  be  the  most  accurate  of  the  three 
methods  studied.  Despite  this,  Truckenbrodt's  formula  was  preferred  over 
the  other  two  methods  due  in  large  part  to  its  simplicity.  It  seemed 
unlikely  to  the  author  that  any  overall  conclusions  drawn  on  the  basis  of 
this  formula  would  be  drastically  overturned  by  use  of  the  momentum 
integral  method.  The  experimental  data  were  obtained^Gertler  ^for  a 
series  of  eight  bodies  of  revolution  in  a towing  tank.  Results  of  the 
comparison  are  shown  in  Table  4.  The  comparison  is  shown  in  terms  of  the 
r.m.s.  fluctuations  between  theory  and  experiment.  No  results  are  reported 
for  the  momentum  integral  approach  at  a Reynolds  number  of  ten  million. 


Nr  = 

Np  = 2 X 10^ 

Truckenbrodt 

3% 

5% 

finite  difference 

6% 

8% 

momentum  integral 

— 

3% 

TABLE  4 - RMS  % ERROR  IN  HESS'  DRAG  COMPUTATIONS  (See  Ref.  27) 


The  small  number  of  bodies  used  by  Hess  to  justify  the  use  of  the 
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Tuckenbrodt  formula  was  deemed  insufficient  by  White  to  establish  the 
method's  validity.  An  extensive  comparison  of  the  Truckenbrodt  formula 
with  experimental  drag  data  obtained  from  seven  series  of  model  forms 
comprising  nearly  fifty  bodies  was  undertaken.  The  computed  drags  show 
very  little  sensitivity  to  such  parameters  as  the  prismatic  coefficient, 
the  position  of  maximum  section,  and  the  nose  and  tail  radii.  The  author 
concludes  that  the  Truckenbrodt  formula  cannot  be  recommended  for 
preliminary  design  or  drag  evaluation  purposes,  and  that  while  the  formula 
does  in  certain  cases  predict  drag  trends  which  agree  with  experimental 
data,  it  is  by  no  means  capable  of  defining  optimums  and  is  very  likely  to 
be  in  error  for  many  cases.  A summary  of  White's  comparisons  appears  in 
Figure  28. 
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Kerney  and  White  recently  introduced  a method  for  the  calculation 

of  viscous  drag  on  streamlined  (i.e.,  non  separating)  bodies  of  revolution. 

Their  approach  was  to  use  an  integral  method  based  on  Granville's, 

supplemented  by  recently  improved  methods  for  predicting  the  transition 

point  and  for  calculating  the  turbulent  boundary  layer.  The 

Douglas- Neumann  program  is  used  to  compute  the  inviscid  pressure 

distribution  over  the  body.  This  is  used  as  an  input  to  the  boundary  layer 

calculations  which,  in  the  laminar  case,  are  identical  to  those  of 
30 

Granville.  The  transition  point  is  found  from  a new  correlation  of 

transition  data  for  bodies  of  revolution  in  terms  of  the  difference  between 

the  momentum- thickness  Reynolds  number  at  the  transition  and  neutral 

stability  points,  and  the  rate  of  change  of  body  shape.  This  method  was 

34 

proposed  by  Granville  and  applied  to  curved  bodies.  Their  predictions 
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son  for  Computed  and  Experimental  Residuary  Drag  Coefficients 


were  compared  with  data  on  bodies  of  revolution  with  and  without  parallel 

middlebodies  from  recent  experiments.  Bodies  without  straight  middle 

sections  were  represented  by  the  DOLPHIN,  from  drop  tests  in  the  Pacific 
35  36 

Ocean,  * whereas  data  on  three  bodies  of  revolution  with  parallel 

middlebodies  came  from  recent  DTNSRDC  towing  tank  tests  of  models  4620-2, 

4620-3  and  4620-4.  Predicted  and  observed  drag  coefficients  and  transition 

points  for  DOLPHIN  exhibited  considerable  scatter,  which  was  attributed  to 

variations  in  the  ocean  ambient  turbulence  level  between  tests,  but 

Granville's  1953  method  was  deemed  a better  predictor  of  transition. 

Consequently,  drag  was  also  predicted  more  acc  iitely  with  this  method. 

Models  4620-2,  4620-3,  and  4620-4  have  streamlined  forebodies  and 

afterbodies  separated  by  parallel  middlebodies  (see  Figure  29).  Drag 

predictions  in  the  3%  to  1%  accuracy  range  were  typical,  with  better 

accuracies  obtained  above  Reynolds  numbers  greater  than  ten  million.  The 
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authors  conclude  that  the  updated  turbulent  boundary  layer  of  Granville 

is  capable  of  giving  reasonable  predictions  of  the  viscous  drag  on 

streamlined  bodies  of  revolution  as  long  as  the  transition  point  is 

accurately  located.  Their  program,  however,  does  not  take  into  account 

displacement  effects  (which  would  increase  the  form  drag).  For  Reynolds 

numbers  less  than  ten  million,  they  suggest  that  a better  method  of 

prediction  for  the  transition  point  is  needed. 
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A.M.O.  Smith  and  T.  Cebeci  presented  the  results  of  an  attempt  to 
improve  presently  available  methods  for  predicting  viscous  drag.  Two 
approaches  were  taken.  In  the  first,  they  attempted  to  directly  compute 
the  drag  from  the  shear  and  pressure  stresses  with  the  Douglas-Neumann 


potential  flow  program,  and  the  Cebecl-Smith  turbulent  boundary  layer 
program,  which  In  principle  are  Ideally  suited  for  this  type  of  study. 

This  approach  was  found  to  produce  reasonable  shear  stresses,  but  the 
pressure  drag  calculations  were  all  wrong,  being  negative  in  some 
instances.  The  exact  reason  for  this  strange  result  is  not  known,  but  the 
authors  suspected  difficulties  arising  partly  from  improper  boundary 
conditions  near  the  body's  trailing  edge. 

The  second  approach  was  to  stick  to  the  basic  momentum  defect  method, 
and  to  try  to  improve  the  calculations  in  the  wake,  resulting  in  a possible 
improvement  over  the  Squire-Young  formula.  Direct  numerical  calculations 
(finite  difference)  were  coupled  to  the  Squire-Young  formula.  The 
Squire-Young  formula  was  used,  but  the  boundary  layer  parameters  were 
obtained  from  finite  difference  calculations  using  the  Cebeci-Smith  method. 
The  turbulent  calculations  on  the  body  were  done  using  the  eddy  viscosity 
concept,  and  the  calculations  were  extended  into  the  wake  by  using  the 
mixing  length  approach.  This  enabled  the  computation  to  proceed  through 
the  problematic  region  of  large  changes  near  the  trailing  edge.  The 
results  were  found  to  be  closer  to  experimental  data  than  the  unmodified 
Squire-Young  approach,  leading  the  authors  to  conclude  that  their  method  is 
promising.  They  noted,  however,  large  differences  in  their  results 
depending  on  how  far  downstream  the  wake  calculation  was  carried.  Much 
more  work,  therefore,  needs  to  be  done  in  this  area  before  a design 

technique  can  take  advantage  of  the  modified  Squire-Young  technique. 
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Parsons  et  al . presented  a method  for  the  synthesis  of  minimum  drag 
hulls  for  axisymmetric  vehicles  of  specified  enclosed  volume  and  constant 
speed  submerged  In  Incompressible  viscous  fluid  at  zero  incidence.  The 
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flow  was  assumed  to  be  nonseparating  and  noncavitatlng.  Drag  reduction  was 
accomplished  solely  through  variation  of  vehicle  hull  shape.  They  claimed 
their  drag  model  to  be  representative  of  the  state-of-the-art  in  drag 
prediction  methods. 

The  boundary  layer  was  assumed  laminar  from  the  stagnation  point 
(leading  edge)  up  to  a transition  point  after  which  it  was  assumed  to  be 
fully  turbulent.  Turbulent  separation  was  not  allowed.  The  mainstream 
inviscid  velocity  distribution  along  the  body  surface  was  computed  with  the 
help  of  the  Douglas- Neumann  program;  viscous  displacement  effects  were 

not  considered.  The  incompressible  axisymmetric  boundary  layer  equations 
were  solved  iteratively  using  an  implicit  finite  difference  method 
developed  by  Cebeci  and  Smith.  A two- layer  eddy  viscosity  model  is  used 
in  the  turbulent  boundary  layer  computation.  Boundary  layer  transition  to 
turbulence  is  predicted  by  the  Smith-Gamberoni  method,  or  by  laminar 
separation/assumed  turbulent  reattachment,  whichever  occurs  first.  The 
Smith-Gamberoni  method  calculates  from  the  Orr-Sommerfeld  equation  the 
amplification  rate  of  Tollmien-Schlichting  waves  in  the  laminar  boundary 
layer.  Transition  is  assumed  to  occur  where  the  amplification  factor 

9 

exceeds  a certain  empirically  determined  value,  in  this  case.e  . 

The  drag  coefficient  is  calculated  using  Young's  formula^® 

„ (H  + ^/2 
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evaluated  at  the  trailing  edge,  where  V is  the  body  volume,  and  u^  is  the 
velocity  at  the  edge  of  the  boundary  layer. 


No  direct  comparison  with  experimental  data  is  presented,  however,  a 
design  is  arrived  at  which  has  a lower  drag  coefficient  (0.0051)  than  the 

DOLPHIN  (0.0076)  at  Re  = 10^. 
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In  1970,  W.L.  Moore  reported  a computer  program  bringing  together 
the  best  techniques  then  available  for  calculating  the  forces  on  the 
components  of  a hydrofoil  boat.  All  flows  were  considered  turbulent,  with 
the  Schoenherr  mean  line  used  throughout  the  calculations  to  define  the 
friction  coefficient,  since  a compromise  had  to  be  effected  between 
accuracy  and  required  computing  time.  The  techniques  used  to  calculate 
strut  drag  are  those  suggested  by  Hoerner.^^  The  foil  profile  drag  is 
calculated  by  using  the  formula 

Cj  - 2 Cf  [l  t 2 (i)  <■  60  (^)  J 


where  is  the  Schoenherr  friction  coefficient,  and  t/c  is  the  maximum 
thickness-to-chord  ratio  on  the  foil.  This  formula  includes  the  effects  of 
skin  friction  and  form  drag  on  the  foil.  This  formula  is  used 
incrementally  and  then  integrated  numerically  to  account  for  any  variations 
of  strut  thickness  ratio  and  strut  chord  along  the  strut  length.  In  the 
case  of  surface-piercing  struts,  an  empirical  term  is  added  to  account  for 
spray  drag,  with  the  spray  drag  coefficient  equal  to  0.24.  For  totally 
submerged  struts,  a term  is  included  to  account  for  the  interference  drag 
between  strut  and  hull 

t t 2 

Cp  = 0.75  (|)  - 0.0003  (|) 

where  the  reference  area  is  taken  to  be  the  square  of  the  strut  thickness. 


85 


itiihll 


Pod  drag  is  calculated  using  empirical  equations  fitted  to  certain 
unpublished  TMB-series  58  body  of  revolution  drag  measured  at  various 
angles  of  attack.  No  further  details  are  given  on  pod  drag. 

Strut-pod  interference  drag  is  calculated  using  Hoerner's  method, 
which  is  to  calculate  twice  the  drag  of  that  portion  of  the  strut  covered 
by  the  pod.  This  is  added  to  the  calculated  drag  of  the  uncovered  part  of 
the  strut. 

The  method  was  primarily  developed  for  strut-foil  and  strut-pod-foil 
combinations,  as  reflected  in  the  awkward  presentation  of  the  drag  data: 
the  lift  coefficient  is  plotted  against  the  drag  coefficient.  Drag 
predictions  are  good  at  low  values  of  the  lift  coefficient  but  rapidly  drop 
b«l  ow  the  10%  accuracy  level  at  lift  coefficients  corresponding  to  angles 

of  attack  exceeding  6 degrees. 
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Nakayama  and  Patel  calculated  the  viscous  resistance  of  bodies  of 
revolution  in  the  absence  of  separation.  Their  method  consists  in 
separately  calculating  laminar  and  turbulent  components  of  the  resistance, 
with  the  transition  point  to  be  determined  empirically.  Although  their 
paper  emphasizes  difficulties  associated  with  the  thickening  of  the 
boundary  layer  near  the  tail  of  the  body,  their  drag  calculations  are 
relevant  to  our  study.  They  find  that  the  resistance  of  bodies  of 
revolution  can  be  predicted  with  considerable  accuracy  provided  the  exact 
location  of  the  transition  point  is  known  or  can  be  found.  This  gives 
confidence  in  the  various  integral  methods  which  are  used  to  calculate  the 
viscous  resistance  in  various  flow  regimes  (laminar,  turbulent),  and 
emphasizes  the  importance  of  being  able  to  predict  the  transition  from  one 
flow  regime  to  another. 
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From  a potential  flow  calculation,  Nakayama  and  Patel  proceed  to 
compute  the  boundary  layer  momentum  thickness  with  a simple  quadrature 
formula  of  the  form 
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which  Is  essentially  the  method  of  Thwaltes.  Here,  6.  is  the  momentum 

thickness,  v the  kinematic  viscosity,  x the  distance  measured  along  the 

surface  from  the  stagnation  point,  r^(x)  the  radius  distribution  of  the 

body,  and  Ug(x)  Is  the  freestream  (potential)  velocity  outside  the  boundary 

layer.  The  momentum  thickness  was  used  in  the  prediction  of  the  transition 

to  turbulence.  Other  boundary  layer  parameters  required  in  the  transition 

prediction  were  obtained  from  Pohl hausen-type  velocity  profiles. 

Transition  to  turbulence  was  determined  from  four  empirical  or 

44  30 

semi -empirical  criteria,  namely  those  of  Michel,  , Granville, 

Crabtree,^®  and  van  Driest  and  Blumer.^®  The  turbulent  boundary  layer 

29 

calcuU'tlons  were  carried  out  according  to  a method  developed  by  Patel. 

The  drag  was  evaluated  by  a computation  of  the  momentum  deficit  at  the  tall 
of  the  body,  by  using  Young's  formula. 

Comparisons  with  experimental  data  revealed  that  no  one  transition 
prediction  method  Is  to  be  preferred  over  the  other.  In  view  of  this, 
experimentally  obtained  transition  point  locations  were  used  In  boundary 
layer  calculations.  In  other  words,  rather  than  using  empirically  derived 
transition  Information,  the  authors  sought  guidance  from  experiments  for  i 

i 

this  part  of  the  computations.  It  was  subsequently  found  that  the  | 


1 


development  of  the  laminar  and  turbulent  boundary  layers  Is  predicted 

satisfactorily,  provided  the  transition  point  Is  accurately  specified. 

Resulting  drag  predictions  were  found  to  be  quite  good,  many  within  one 

percent  of  the  experimental  value. 
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Cebeci  et  al.  presented  calculations  of  viscous  drag  for  two 
dimensional  and  axisymmetric  bodies.  Drag  was  calculated  from  Granville's 
formula 


r _ 20 
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evaluated  at  the  trailing  edge  of  the  body.  The  above  formula  reduces  to 

the  Squire-Young  formula  for  two  dimensional  bodies  by  setting  n » 1 In 

which  case  e Is  the  momentum  thickness  and  A Is  the  chord.  For 

axisymmetric  bodies,  n * 7,  © is  the  momentum  area,  and  A is  a suitable 

area  characteristic  of  the  body.  Cebeci  et  al.  obtained  their  boundary 

layer  parameters  not  from  an  Integral  approach,  but  from  a more  exact 

22 

finite  difference  solution  developed  by  Cebeci  and  Smith.  Boundary  layer 

separation  was  not  considered.  To  predict  transition  to  turbulence,  the 

44 

empirical  method  of  Michel  was  utilized.  This  method  is  based  on  a 
correlation  of  transition-momentum-thickness  Reynolds  number  ^ with 
transition  x-Reynolds  number,  Np  approximated  by  the  formula 
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I I Theoretical  predictions  of  the  drag  coefficients  were  compared  to 

I experimental  data  for  airfoils.  The  r.m.s.  error  based  on  57  drag  values 

[ ! ’ was  found  to  be  2.9X.  The  error  In  calculated  drag  Increased  with  the 

I i 

angle  of  attack.  Displacement  thickness  corrections  were  found  necessary 
for  bodies  with  fineness  ratios  {length/maximum  diameter)  less  than  4, 

I while  for  slenderer  axlsymmetric  shapes,  skin  friction  accounted  for  the 

[ experimental  drag. 

DRAG  PREDICTION  METHODS  ACTUALLY  USED  IN  DESIGN  PROGRAMS 

In  this  section  we  discuss  drag  prediction  methods  which  have  been 
Implemented  In  strut  pod  design  programs. 

Developmental  Sciences,  Inc.  (DSI)  designed  a pod  1nlet^  for  a 200 
ton,  100  kt  hydrofoil  craft.  The  external  pod  shape  was  determined  by  the 
requirement  that  there  be  no  external  cavitation  at  cruise  speed.  Details 
about  their  calculation  methods  are  sparse,  making  an  evaluation  of  their 
design  difficult.  The  external  drag  was  determined  by  computing  the 
turbulent  boundary  layer  over  the  exterior  surface  and  Integrating  the  skin 
friction  up  to  the  point  of  separation,  the  calculations  were  then 
continued  beyond  the  separation  point  by  using  an  unexplained  empirical 
relationship.  Four  options  were  analyzed:  1)  fully  wetted  pod  with  short 
boattall,  11)  a linear  extension  of  the  boattall  which  eliminated 
separation.  111)  a ventilated  base  and  1v)  extended  base  ventilation.  The 
relative  magnitudes  of  skin  friction  and  form  drag  are  not  discussed,  nor 
Is  any  account  made  of  the  strut  drag,  or  strut-pod  Interference  effects. 


1 


89 


A very  understandable  account  of  the  problems  associated  with 

3 

strut-pod  Inlet  design  was  presented  by  Sherman  and  Lincoln.  The 
stagnation  streamline  divides  Internal  and  external  flows  and  the  total 
drag  Is  the  sum  of  Internal  and  external  drags. 


^NET  ' *'lNT  '"ext* 
The  external  drag  coefficient  Is  given  by 
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where  Is  the  Inlet  area,  V^/V^  the  IVR,  and  1s  the  "system  drag" 

which  Is  the  summation  of  all  the  external  viscous  drag  forces  plus  any 
external  pressure  and  Interference  drags  not  previously  accounted  for.  The 
second  term  In  the  above  expression,  dubbed  "pre-entry  drag,"  has  its 
origin  In  the  acceleration  of  the  flow  In  the  region  ahead  of  the  Inlet. 

For  fully  settled  turbulent  flow,  the  system  drag  coefficient  of  the 
ram  Inlet  nacelle  based  on  Its  external  wetted  area  Is  estimated  from  the 

following  formula  allegedly  obtained  from  Hoerner;^^ 

[d„1.5  d^  S'! 

where  d^^^  Is  the  maximum  external  nacelle  diameter,  the  nacelle  length 
and  Is  not  defined  by  the  authors.  Similarly,  the  fully  wetted 
external  turbulent  viscous  drag  coefficient  for  the  supporting  strut  Is 
obtained  from  the  same  source: 

^0  ' 2 'fs  [’  * 2 ] 
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where  Is  the  thickness  to  chord  ratio.  Any  Interference  drag  between 
the  nacelle  and  strut  Is  considered  negligible  for  small  values  of  strut 
thickness  to  chord  ratio.  The  spray  drag  based  on  strut  planform  area  Is 
given  by 


for  Froude  numbers  (based  on  the  chord)  greater  than  3. 

The  external  viscous  drag  of  the  system  Is  reduced  If  a base  vented 
strut  and  nacelle  are  used.  In  this  case  spray  drag  may  be  anywhere  from 
1/6  to  1/2  of  the  total  strut  drag.  Using  Hoerner's  method  for  determining 
the  total  drag  of  a base  vented  strut  (Including  spray  drag),  the  following 
formula  Is  obtained 


Again,  the  nomenclature  Is  not  explained.  The  total  drag  of  a base  vented 
nacelle  (based  on  the  wetted  nacelle  area)  Is  given  as 


^ 9 "n  ' 'm 


where  g Is  the  acceleration  due  to  gravity,  h^  Is  the  submergence,  V Is  the 

speed,  A|^  the  wetted  area  and  the  maximum  radius  («  d^^/Z). 

18 

The  strut  drag  In  supercavltating  flow  Is  given  as 
Cjj  (a)  » Cp  (o  - 0)  (1  + o) 
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where 


P - P 


cav 
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Is  the  cavitation  number,  and  Cp{o)  is  the  drag  coefficient  at  zero 
cavitation  number.  The  authors  produce  the  following  formula  for  the 
system  drag  coefficient  for  a supercavltating  ram  Inlet 


■ “N  (O'  (’  * 


where  the  first  term  Is  the  nacelle  cavitation  number,  and  the  second  term 
represents  strut  drag  based  on  inlet  area. 

For  the  purpose  of  internal  loss  analysis  of  waterjet  inlets  the 
inlet  Is  broken  into  three  components:  the  inlet  diffuser,  the  transition 
elbow  Joining  the  strut  and  nacelle,  and  the  strut  diffuser  section.  A 
total  head  loss  is  obtained  for  the  inlet  by  summing  the  component  losses. 
Elevation  losses  are  not  considered  in  these  calculations. 

The  inlet  diffuser  is  a critical  component  of  the  waterjet  duct 
system.  Since  the  dynamic  pressure  is  greatest  at  the  diffuser  entrance, 
an  inefficient  diffuser  will  seriously  degrade  overall  propulsor 
performance  at  high  speeds.  Sherman  and  Lincoln's  method  of  predicting 
diffuser  losses  Is  an  approximate  two-dimensional  laminar  and  turbulent 
boundary  layer  calculation  beginning  at  the  inlet  stagnation  point  and 
terminating  at  the  maximum  Internal  inlet  diameter.  The  effects  of  axial 
syimietry  and  wall  curvature  on  boundary  layer  development  are  not 
considered. 

The  Neumann  problem  Is  solved  to  obtain  the  potential  flow.  A 
step-by-step  solution  of  the  Integral  boundary  layer  equations  is  then  used 
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to  obtain  the  distribution  of  momentum  thickness .displacement  thickness  and 
shape  factor  along  the  diffuser. 

Transition  Is  assumed  to  occur  when  a certain  critical  value  of  the 

Op 

displacement  thickness  Is  reached.  The  turbulent  portion  of  the  boundary 
layer  Is  calculated  from  Truckenbrodt's  Integral  analysis.  Starting  with 
the  value  of  momentum  thickness  at  the  transition  point,  the  momentum 
thickness  turbulent  shape  factor,  and  the  displacement  thickness  are 
calculated  along  the  wall.  The  value  of  the  shape  factor  for  which 
separation  will  occur  Is  generally  found  to  be  between  1.8  and  2.4. 

However,  a shape  factor  of  3.1  has  been  shown  applicable  to  conical 
diffusers  and  Is  used  as  the  upper  limit  1n  Sherman  and  Lincoln's  analysis. 

With  subscripts  1 and  2 denoting  the  diffuser  entrance  and  exit, 
respectively,  the  diffuser  pressure  recovery  efficiency  Is  calculated  as 
follows.  Once  the  displacement  thickness  Is  known,  the  effective  area  at 
the  diffuser  exit  Is  given  by 

^?^eff  ' (^2  ■ 

The  overall  static  pressure  recovery  Is  defined  as 


P2  - Pi 

q 


q = >5  p V 


1 


The  Ideal  static  pressure  recovery  for  Inviscid  one  dimensional  flow  Is  | 

given  by  1 


'pr.l 


A,  2 
^2 


93 


The  diffuser  efficiency  or  effectiveness  Is  simply  the  ratio  of  overall  to 
Ideal  pressure  recovery. 


pr 


.1 

q 1 - 


^2 


Assuming  a one  dimensional  potential  core  flow  outside  the  boundary  layer, 

the  diffuser  efficiency  becomes 

A,  2 
1 - (^/ 

'^2  eff 


nn  = 


A,  2 
1 - (^) 


The  total  head  loss  In  the  nacelle  elbow  Is  expressed  by  an  equation 
of  the  form 


e ^ 


where  the  friction  factor  f Is  obtained  from  circular  pipe  data  using  the 
equivalent  hydraulic  diameter  d^,  and  L Is  the  elbow  length.  With  complex 
geometries,  component  losses  are  assumed  to  be  additive.  The  effect  of  channel 
curvature  on  the  friction  factor  Is,  in  the  above  equation,  given  by  replacing 
f by 


where  r Is  the  turning  radius.  The  straight  duct  friction  factor  f Is 
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calculated  from  Nikuradse's  formula 


^ = 4 log^Q  (Np  /r  ) - 0.40 

The  strut  channel  Is  treated  as  a rectangular  diffuser  with  two 
convergent  walls.  Again  using  the  upstream  and  downstream  subscripts  1 and 
2,  respectively,  strut  losses  are  calculated  using  the  following 
approximation  for  a truncated  diffuser: 


'l=  [' '4’] 


where  the  friction  factors  are  given  In  Reference  34. 

The  Lockheed  design  method  also  resorts  to  a formulary  for  the 
calculation  of  drag.  In  the  subcavltating  case,  strut  drag  Is  composed  of 
parasite  and  spray  drag.  Strut  parasite  drag  1s  estimated  by 

Cp  = 2 Cf  1^1  + 1.2  (^)  + 60 

where  Is  the  Schoenherr  friclton  factor.  Spray  drag  Is  obtained  from 
Cp  = 0.24  t^ 

2 

where  = >sVq  In  the  dynamic  pressure. 

Without  further  ceremony,  the  Lockheed  report  refers  the  reader  to 
Hoerner^^  and  "various  NACA  publications"  for  method  of  estimating  strut, 
foil  and  nacelle  Interference  drag.  The  drag  coefficient  for  the  nacelle 
Is  estimated  by  the  following  equation  which  Is  based  on  wetted  area 

3/2  d„  31 
1 + 1.5  (j^)  + 7 (j^)  J 


where  d_  Is  the  maximum  external  nacelle  diameter  and  L Is  the  maximum 
n 

nacelle  length. 

The  drag  of  base-vented  nacelles  Is  determined  by  the  differences  In 
pressure  distribution  between  the  base  and  the  nose  after  deducting  the 
Internal  momentum  drag.  The  nose  pressure  Is  the  sum  of  atmospheric  and 
hydrostatic  pressures.  While  the  base  pressure  Is  atmospheric.  The  total 
base  drag  Is  then  taken  to  be  the  product  of  hydrostatic  pressure  times  the 
base  area. 

Interference  drag  of  strut,  nacelle  and  foil  Intersections  resulting 
from  flow  separation  along  the  afterbody  Intersections  Is  considered 
negligible. 

Internal  pressure  losses  are  caused  by  a)  friction  on  the  duct  walls 
and  b)  form  effects  due  to  changes  in  the  size,  shape  or  direction  of  the 
duct  and  are.  In  general,  proportional  to  the  kinetic  energy  of  the  fluid 
for  any  duct  configuration.  For  circular  duct,  the  surface  area  Is  trLd 
(L  = duct  length,  d = diameter)  and  the  pressure  loss  can  be  expressed  In 
terms  of  the  Fanning  friction  factor  F 


where  q Is  the  dynamic  pressure,  d^,  the  hydraulic  diameter,  is  defined  by 

d a 4 cross  sectional  area 
h ^ wetted  perimeter 

and  F Is  given  graphically  In  the  Lockheed  report  as  a function  of  Reynolds 
number  and  surface  roughness.  For  all  ducts  associated  with  waterjet 
systems,  Lockheed  recommends  that  the  smooth  pipe  curve  be  used. 
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I 


Pressure  losses  due  to  changes  In  shape  and  direction  of  the  duct  are 


proportional  to  the  dynamic  pressure  and  can  be  expressed  as 

where  ICj.  depends  on  the  particular  duct  configuration.  The  total  pressure 
loss  for  any  component  can  be  expressed  by  the  sum 


q 


K.  + 4 F V 
^ '‘h 


The  purpose  of  diffuser  Is  to  convert  kinetic  energy  at  the  entrance 
Into  pressure  at  the  exit  with  a minimum  loss  In  total  pressure.  The  total 
pressure  loss  Is  directly  related  to  the  Internal  drag  of  the  component 
while  the  exit  flow  distribution  affects  the  operation  of  the  duct 
components  located  downstream  of  the  diffuser.  It  has  been  shown  (although 
no  reference  is  given  in  the  Lockheed  report)  that  for  a given  expansion 
angle  and  fixed  Inlet  conditions,  the  total  pressure  loss  of  a diffuser  Is 
proportional  to  the  theoretical  value  of  total  pressure  loss  for  a sudden 
expansion  of  the  same  area  ratio.  The  proportionality  factor  is  defined  as 


where  Ap  = A2/A^.  and  the  subscripts  1 and  2 define  respectively  the 
diffuser  Inlet  and  exit  stations,  denotes  the  total  pressure  and  q the 
dynamic  pressure. 
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Two  geometric  variables  are  considered  Important  In  diffuser 
operation:  area  ratio  and  expansion  angle.  In  general,  Lockheed  considers 
that  the  design  of  a good  waterflow  diffuser  can  be  accomplished  by  having 
a uniform  Inlet  velocity  profile  and  keeping  the  expansion  angle  2a  between 
8 and  12  degrees.  It  Is  also  advantageous  to  prevent  local  flow  separation 
due  to  wall  roughness  and  sudden  changes  In  flow  direction.  If  possible, 
the  diffuser  axis  should  be  straight.  It  Is  also  claimed  (again,  no 
reference  cited)  that  diffusers  having  an  area  variation  such  that  the 
axial  pressure  change  varies  directly  with  the  length  achieves  a better 
exit  total  pressure  recovery  than  a conical  diffuser  of  the  same  length  and 
contraction  ratio.  A diffuser  of  this  type  has  a so-called  "trumpet"  shape 
with  an  increasing  rate  of  change  of  cross  sectional  area  with  distance 
from  the  Inlet.  (Interestingly,  the  trumpet  shape  is  condemned  by  Sherman 
and  Lincoln  as  particularly  susceptible  to  separation). 


= constant 


The  loss  coefficient  for  bends  and  elbows  is  a function  of  many 
variables:  duct  cross  sectional  shape,  turning  angle,  turning  radius,  area 
variation  and  Reynolds  number.  The  total  loss  coefficient  for  elbows  is 
given  by 


^ Pt 
q 


C K.  + 4 F L 
*90  3;; 


where  K^.  is  the  loss  factor  for  90  degree  bends  and  c is  a correction 
^90 

factor  for  bends  other  than  right  angle.  These  quantities  are  presented 


graphically  In  the  Lockheed  report.  Losses  In  ducts  may  be  reduced 
considerably  by  employing  splitters  or  turning  vanes  which  divide  the  flow 
In  parallel  channels,  each  having  a larger  radius  to  diameter  ratio  and 
aspect  ratio  than  the  original  elbow.  Care  must  be  exercised  to  Insure 
that  the  Increased  friction  losses  resulting  from  the  splitters  are  not 
greater  than  the  Improvement  In  form  loss. 


REVIEW  OF  THE  EXPERIMENTAL  DATA 

Presently  attainable  drag  performance  levels  are  Indicated  In  Table 
2.  The  drag  coefficient  (based  on  Inlet  area  and  cruise  speed)  of  the 
Aerojet  Inlet  Is  0.45  and  not  Inordinately  out  of  line  with  the 
representative  figure  of  0.35  given  by  EtterJ^  Sobolewski's  measurements® 
give  1.67  for  the  DSI  model  Inlet  system,  but  It  must  be  remembered  that 
this  figure  Includes  strut  drag. 

The  predictive  accuracy  of  existing  drag  prediction  methods  Is 
described  In  Table  5. 

In  studying  waterjet  Inlet  technology,  one  becomes  acutely  aware  of 
the  deficiencies  In  a field  which  has  not  yet  reached  maturity.  In  the 
present  situation,  one  faces  a variety  of  experimental  measurements  of 
"Inlet  drag"  by  Investigators  who  do  not  subscribe  to  the  same  definition 
of  this  Important  quantity.  Although  all  drag  measurements  reported  here 
were  made  on  ram  Inlets,  some  data  include  foil  drag  and  others  contain  the 
effects  of  Internal  ducting  and  even  the  exit  nozzle.  Strictly  speaking. 
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Method 

Investigator 

% Error 

, ! 

Truckenbrodt 

Hess  (27) 

4 

♦ 

Momentum  Integral 

Hess  (27) 

3 

. j 

Young's  Formula 

Hess  (27) 

7 

1 

i 

i 

Truckenbrodt 

White  (32) 

>50  (see  text) 

! 

Granville 

Kerney,  White  (33) 

3-7 

1 

Hoerner 

Moore  (41) 

*10 

Nakayama-Patel 

Nakayama,  Patel  (42) 

1 

Cebeci 

Cebeci  et  al.  (47) 

3 

TABLE  5 - APPROXIMATE  PREDICTIVE  ACCURACY  OF  DRAG  PREDICTION 
METHODS 

each  measurement  Is  therefore  valid  only  for  the  particular  hardware  on 
which  It  was  performed.  Fundamental  differences  between  drag  measurements 
performed  on  various  Inlet  concepts  preclude  a direct  comparison  of  their 
relative  merits.  Furthermore,  one  must  regard  the  available  data  with  a 
skeptical  mind. 

The  effect  of  the  variable  area  factor  (the  ratio  of 
maximum  to  minimum  Inlet  area)  on  Inlet  drag  has  been  presented  by  Etter.^^ 

Smaller  Inlet  drags  are  obtained  with  large  variable  area  factors  since  ^ 

less  leading  edge  (lip)  thickness  Is  needed.  In  other  words,  by 
accommodating  a demanding  flowrate  schedule  with  a variable  area  Inlet,  one 
no  longer  needs  large  and  draggy  lip  thickness.  The  trade-off  between  lip 


I 

t 


thickness  and  variable  area  factor  is  such  that  an  overall  lesser  drag  is 

't 

obtained  with  variable  area  inlets.  While  Etter's  conclusions  are 
demonstrated  for  flush  inlets,  there  is  no  reason  to  expect  a different 
conclusion  for  strut  pod  inlets. 

CONCLUSION  AND  RECOMMENDATION  CONCERNING  INLET  DRAG 

^ In  reviewing  drag  prediction  methods,  it  became  clear  that  particular 

attention  must  be  paid  to  the  various  flow  regimes  (laminar,  turbulent, 
j separated)  in  order  to  arrive  at  an  accurate  drag  figure.  Various 

' approximate  approaches  have  been  developed  to  identify  these  regimes  and  to 

compute  the  resulting  drag.  The  method  of  Nakayama  and  Patel  appears  to 

be  quite  accurate  (1%  error)  and  utilizes  only  flow  integral  parameters, 

! * meaning  that  a detailed  knowledge  of  the  boundary  layer  structure  is  not 

needed.  It  must  be  said  that  this  is  an  active  field  of  research  and  that 
more  literature  on  this  subject  is  likely  to  be  generated  in  the  near 
future,  but  a drag  prediction  technique  based  on  an  integral  flow  analysis 
appears  to  be  promising  for  inlet  applications. 

3 

Sherman  and  Lincoln  have  demonstrated  the  importance  of  internal 
I losses  in  the  design  of  strut-pod  inlets.  A reliable  method  of  estimating 

these  losses  in  strut  pod  structures  is  needed,  and  can  be  formulated  from 
an  empirical  study  of  pressure  profiles  in  such  structures. 
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RECOffCNDATIONS  FOR  FURTHER  WORK 


In  view  of  the  eventual  goal  of  establishing  a reliable  design  method 
for  strut-pod  Inlets,  special  efforts  should  be  concentrated  In  the 
following  areas. 


[ 1.  It  should  be  determined  If  the  nacelle  and  Inlet  shapes  proposed  by 

i 

I Sherman  and  Lincoln  are  Indeed  optimal  from  a cavitation  point  of  view. 

[ Systematic  ways  of  producing  cavitation  free  nacelle  contours  should  be 

[‘ 

i ) Investigated. 


2.  There  Is  an  urgent  need  for  a definitive  report  on  Inlet  drag  and 
losses.  This  would  Include  analytical  work  aimed  at  clearly  Identifying 
the  various  components  of  drag  and  Internal  losses,  their  effect  on  Inlet 
and  craft  performance,  and  would  serve  as  a guide  to  experimenters  In 
making  measurements  that  can  be  meaningfully  compared. 

3.  As  previously  mentioned,  the  Sherman-LIncoln  design  philosophy  of 
minimizing  overall  losses  In  cavitation- free  Inlets  should  be  coded  as  a 
modular  program. 

4.  A review  of  optimization  techniques  should  be  made  to  find  a rapid  and 
efficient  method  of  evaluating  the  many  trade-offs  encountered  In  Inlet 
design.  Specifically,  a numerical  optimization  scheme  should  be 
recommended  for  Inclusion  In  the  modular  design  program. 
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5.  A review  of  presently  existing  pressure  data  obtained  from  tests 
performed  on  actual  Inlets  should  be  examined  with  the  objective  of 
formulating  a reliable  Internal  loss  formula  based  on  actual  model  tests. 
Such  data  already  exist  and  If  needed,  further  tests  could  be  performed  on 
the  (suitably  modified)  DSI  pod,  a model  of  which  already  exists  at 
DTNSRDC.  Such  a formula  would  Inherently  Include  Interaction  effects  and 
would  be  specifically  applicable  to  the  strut-pod  concept. 

6.  Another  area  In  need  of  study  Is  the  vane  system  at  the  strut-pod 
elbow.  This  area  Is  susceptible  to  cavitation  and  high  losses.  Several 
approaches  are  possible,  from  a multi-element  vane  system  to  the 
compartmentalizatlon  of  the  elbow  Into  a few  channels  separated  by  flow 
dividers.  A cavitation  and  loss  analysis  should  be  performed  to  develop  an 
optimum  flow  turning  geometry. 

7.  A simple  momentum  analysis  of  the  flow  In  an  around  the  inlet  does  not 
provide  a satisfactory  estimate  of  the  drag,  since  no  control  volume  can  be 
specified  with  known  velocities  and  pressures  on  all  of  Its  surfaces.  An 
alternative  method  is  needed,  based  on  simple  approximation,  to  yield  the 
pre-entry  drag  and  external  pressure  recovery.  This  could  greatly  help 
generate  preliminary  performance  figures  In  a screening  analysis  of 
candidate  Inlet  shapes. 
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